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ABSTRACT 


A  brief  survey  of  previous  developments  in  the  theory  of  partial  coherence 
is  given  in  Chapter  1  and  the  relations  between  the  complex  degree  of  coherence 
(see  Wolf  1955)  and  the  measures  of  coherence  introduced  by  earlier  authors 
are  given.  Because  of  the  sundry  formulations  of  the  theory  to  be  found  in  the 
literature  a  detailed  review  of  the  formalism  used  here  is  also  given  in  this 
chapter. 

Coherence  theory  is  formulated  in  terms  of  correlation  functions  and 
analytic  signals  aiid  the  developments  of  this  thesis  required  several  new  theorems 
concerning  the  convolution  uid  cross-correlation  of  analytic  signals.  These 
theorems  are  developed  in  Chapter  2.  The  form  of  the  general  Lo.imodular 
analytic  signal  is  also  obtained  in  this  chapter.  It  is  shown  that  when  considered 
as  a  function  of  a  complex  variable  this  function  is  a  meromorphic  function  of 
order  one  with  isolated  poles  in  either  the  upper  or  the  lower  half  plane  only 
and  with  zeros  at  conjugate  points. 

In  Chapter  3  these  theoreme  are  applied  to  the  detailed  analysis  of  the 
limiting  forms  of  the  mutual  coherence  function  for  toth  polychromatic  and 
quasl-monochromatlc  fields.  In  the  rigorous  analysis,  applicable  to  fields  of 
arbitrary  spectral  w^dth,  It  Is  shown  that:  1.)  ati  o  'tical  field  is  coherent  if 
and  only  If  it  is  monochromatic;  2.)  the  mutual  coherence  ^unction  for  a  co¬ 
herent  field  may  be  expressed  as  the  product  of  a  simple  weave  function,  evalua¬ 
ted  at  one  point,  with  Its  complex  conjugate,  evaluated  at  a  second  point,  multi- 
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plying  a  simple  periodic  factor,  e  ;  3.'  an  Incoherent  field  cannot  exist 

in  free  space  though  an  incoheren.  source  may  be  defined  in  a  manner  consisstent 
with  this  result.  The  esser.dal  differences  between  these  theorems  and  the 
cor respoiviing  ones  for  qjasi-monochromr tic  fields  are  discussed. 


In  Chapter  4  the  propagation  of  mutual  coherence  is  studied.  A  new 
derivation  of  the  wave  equations  for  this  quantity  is  given  and  the  equations 
are  solved  for  the  field  created  by  an  arbitrary  plane  source,  i.e.,  an  extended 
partially  coherent  polychromatic  source.  Using  the  results  of  the  previous 
chapter,  the  limiting  forms  of  the  general  solution  are  examined  in  detail  for 
both  polychromatic  and  quasi-monochromatic  illumination.  In  particular  it 
is  shown  that ;  1.)  a  coherent  source  always  gives  rise  to  a  coherent  field;  and 
2.)  an  incoherent  source  always  creates  a  partially  coherent  field.  These  re¬ 
sults  are  shown  to  be  valid  regardless  of  the  spectral  width  of  the  illumination. 
By  examining  the  general  solution,  the  well-known  van  Cittert-Zernike  theorem 
is  found  to  be  an  approximate  form  of  the  incoherent  limit  of  the  quasi-mono¬ 
chromatic  solution. 

The  results  of  the  previous  chapters  are  applied  In  Chapter  5  to  a  fre¬ 
quency  domain  analysis  of  the  optical  Imaging  problem.  A  general  solution  to 
this  problem  for  partially  coher.mt  polychromatic  light  Is  obtained  and  genera- 
llzffd  transfer  functions  are  introduced.  It  is  shown  that,  under  the  quasi- 
monochromatic  approximation,  these  generalized  transfer  functions  reduce  to 
the  familiar  forms  found  in  the  literature.  A  transfer  function  is  introduced 
for  obtaining  the  mutual  spectral  density  of  vhe  image,  a  function  of  tw  points, 
from  the  self  spectral  density  of  the  obteci,  a  function  of  one  point. 
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INTRODUCTION 

It  is  cuatomary  to  describe  interference  as5d  diffraction  phanomeria  in 
terms  of  tJje  superposition  of  strictly  coherent  er  strictly  incoherent  vibra¬ 
tions  Ignoring  completely  the  ptissibillty  ’  intermediate  states,  partial  co¬ 
herence.  That  this  practice  is  a  considerably  restrictive  idealization  ia 
attested  to  by  the  fact  that  experience  iixlicates  that  both  of  these  extremes 
are  unrealizable. 

An  emnple  of  the  inadequacy  of  the  concepts  of  complete  coherence  and 
complete  incoherence  for  the  description  of  physically  interesting  phenomena 
was  known  as  early  as  about  1869.  At  that  time  Verdet  (1869)  detnonstraied 
that  the  light  from  two  pinholes  in  a  screen  Illuminated  by  the  sun  will  inter¬ 
fere  la  Young’s  interference  experiment  if  tti«  stparation  of  the  pinholes  is 
less  than  about  l/3S)mm,  Since  interference  is  customarily  a  property  of 
coherent  oscillations,  tiiis  result  jiuggested  the  idea  of  an  "area  oi  coherence” 
su’TOundlng  a  point  in  an  Ineohercni  field. 

Afiothor  early  exampic  of  th-?  inadequacy  of  these  limiting  concepts  can 
be  found  In  the  work  of  Micljelson  from  about  1090  to  a91K).  The  interferometric 
method,  introduced  by  Mlcbelson  (1890-1921),  erf  measuring  stellar  diameters 
esswitUliy  iswolves  measuring  the  d-.‘gree  of  colierence  of  the  lUumlnation 
prcdcced  by  bie  stars.  Although  this  wt>rk  was  not  interpreted  in  terms  of 
coherence  theory  until  mucv’  later,  the  i|uantl'y,  viaibility,  introduced  by 
Michelsoi:  to  describe  the  quHlir;  of  intccferencc  fringes  b4ui  proved  to  be  ai. 
important  key  to  ardor  standing  the  ccncipt  of  parUal  coherence,  tn  fact,  in 
2iei  ntke's  i‘i938)  fcrmulatior'  of ,  C‘';oreoce  theory,  the  degree  of  coherence 


betWMn  the  vlbraticns  at  two  points  in  an  optical  field  is  defined  as  the  visibility 
of  the  fringes  obtained  by  allowing  them  to  interfere  in  a  suitable  eaperiment. 

From  the  middle  ctf  the  nineteenth  century  until  the  iamt  two  decades  the 
theory  of  partial  coherence  received  but  lltUe  attention.  The  few  papers  that 
did  ai^ar  on  the  subject  are  more  or  less  disconnected  since  rich  in**e8tifatox 
introduced  his  own  appar  ently  different  formulation  otf  the  theory.  The  lack  of 
interest  in  the  subject  during  that  period  may  be  attributed  to  the  fact  that 
applications  in  which  the  theory  is  Important  were  either  unknown  or  inve^vad 
measurements  which  were  not  refined  enough  to  take  account  of  the  d^ree  of 
coherence. 

However,  in  more  recent  years  the  concept  oi  partial  coherence  has  be¬ 
come  important  in  virtually  every  branch  of  physics  which  involves  eiectro- 
magnetic  radiation  regardless  of  the  frequencies  considered.  In  visible  i^ties, 
for  example  ,  coherence  theory  is  tantamount  to  the  uMierctanding  of  Kich  topics 
as  image  formation  and  the  effect  of  iUuminitiOii  on  the  reeolution  in  a  micro* 
scepk  image.  bi  spectroecopy  the  influence  of  the  slit  whidi  on  the  degree  of 
coherence  of  the  iUumination  can  produce  measurabie  effects.  In  ra^o  astronomy 
source  dis  leters  are  measured  by  interferoanetric  techniques,  which  eff^tively 
involves  die  meu.surement  of  the  degrse  of  coherence  cf  ihe  radiaUon.  Similar 
proMem*^  arise  in  these  apfdicatione  of  radar  where  ques'^ione  of  resolution 
and  mappbiS  arc  of  central  iutportanc«. 

Even  at  the  much  longer  wave  lengths  used  in  communication,  corierence 
t.Vory  i^ays  an  important  role.  In  scatter  cr-mmunicatione,  for  example,  Beran 
(195S)  has  recently  studiod  the  propikgation  of  the  ensemble  correlation  in  the 
scatter  flrtd  and  wlQi  the  present  author  has  shown  that  ‘he  reliability 

of  a  spatially  diversified  scatter  system  can  be  computed  from  a  knowledge  of 


th«  d^ree  of  cobersnce  in  fhss  reflected  field. 

In  the  last  few  years  as  these  &pf)Ucatlor®  became  unportaat  and  a« 
>si4au»urii^  techndqoes  became  more  refined,  the  th^ry  of  partial  coherence 
hae  received  increaeij^  attention  in  the  fields  mentiot^  aocve.  In  spite  of 
this  aiteniion  several  aspects  cf  the  tl^eor^'  have  been  misufKlerstood;  and  tbs 
strot^  emphasis  on  application  has  left  many  of  the  fundamental  questlo?* 
unansvered. 

Th^  aim  of  this  thesis  is  the  deduction  of  certaU  of  ihe  ger^ral  math«iflftt~ 
leal  aisd  physical  iir.plicati<M5S  cX  coherence  theory.  For  this  purpose  the  a>oet 
convciiieni  and  rigorous  fornialatlon  of  the  subjei’t  is  that  introduced  by  Wolf 
(1^55);  and  feia  formuUmt  will  i>e  rsed  exciusivciy  in  ibis  thesis.  However, 
befere  review in  detail  the  structure  of  the  theory  to  !)«  usiid  here,  it  will 
prove  useful  to  give  a  brii?f  survey  of  the  previcnis  research  on  the  eubjsct, 
pevtoj®  particular  attention  to  the  defUiiHoas  of  defree  of  coherence  iotroduced 
by  eariisr  authors.  In  Shis  way  some  of  the  advantages  of  the  foraiuiaOoe  used 
here  will  become  apparent, 

Kariy  research  on  partial  cohere,<.ce  is  asaoclaied  with  the  names  of  Verdet 
{1-36?J,  von  Laue  (1§07),  Berek  (lt»5  %,  c,  d)  and  van  CUtert  (l!i.S4  !939).  The 
inTestlgftdona  of  this  subject  in  the  last  two  decades  are  found  primarily  in  the 
work  of  Zernike  (1953).  Hopkins  {19SI,  1^57),  Bianf'-Lapierve  and  Dumontet 
{!855),  and  Wolf  (1555).  Since  ^tl  of  the  earlier  formulations  of  the  theory  are 
special  cases  of  the  on*-  used  here,  it  will  prove  convenient  is  introduce  the 
pr'ncipie  fvncbossn  of  this  general  U»sory  before  discuss* ng  th*  contrtoutlona 
s”if  tha  pcfevlous  writers. 
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ThiJ  study  cA  partial  coherence  in  optica  ia  essentially  the  study  dL  the  com- 
piex  cross-cor  relation  of  the  disturbances  at  two  typical  points,  Pj  and  in 
the  optical  field.  This  function  is  defined  as  the  mutual  coherence  fuocticNC*, 

v/  (t))^  , 

where  Vj(t)  and  V2(t)  are  the  "consples**dlaturbaiice8"  at  the  two  points,  U.o 
sharp  brackets,  ,  denote  time  average,  and  r  is  the  delay. 

complex:  degree  of  coherence,  defined  as  the  norcnaliai^d  mutual 

coherence  funcUcn, 


It  should  be  noted  that  the  degree  of  cohei^e  is  a  funcUon  not  orly  of  the 
coordUmtss  of  Pj  and  Pj  but  alao  of  the  time  delay  o’-  path  difference  coor^nate, 
T,  With  these  defiritioos  In  mind  wc  may  now  direct  oui  attention  tc  the  early 
contributions  to  the  subject  of  partial  cc^remre. 

After  the  work  erf  Verdet  (18<J9),  mentioned  wirlier,  a  cjuaniitative  m.»asure 
for  partial  coherence  was  given  In  a  paper  by  von  Laue  (1907).  la  thui 
in  which  was  discussed  the  thermodvnamicai  aspects  of  diffruchon,  the  qua-’?dti’ 

,  propoi’tional  to  square  of  the  time  averajed  product  of  Uie  oisturbaucee 
at  two  points  in  the  field,  proved  to  bo  of  centra!  importance.  While  this  <pii»ntlt>' 
was  sufficient  to  charmc‘ertae  the  optical  field  for  Int*  probleiijc  di8cusa<e^  by 
von  Laue,  the  formulstlon  ie  too  reetrictive  for  general  jqpf’licati&ii.  fiowe-'er ,  ^■'e 
may  mention  that  the  recent  interesting  experiment  of  Ha."bii?y-Brc>wn  sixi  Twles 

•  l^eclse  definitions  id  *hls  and  the  other  furctiaos  relerreu  to  here  are 
given  in  the  next  secUen. 

••  Tbe  preciee  meaning  of  the  term  complex  disturbance  is  given  inSectica  3. 


(1956  a,  b)  measured  precisely  the  quantity  . 
It  is  shown  in  Appendix  1  of  this  tiiesis  that 


The  next  theorotica)  treatment  of  the  subject  appeared  some  twenty 
years  later  in  the  work  of  Berek  (1926)  in  which  was  Introduced  the  so-called 
degree  of  consonance  as  a  measure  of  coherence.  Berek's  formulation  was 
applied  to  some  problems  in  the  theory  of  image  formation  in  the  microscope. 

Some  of  his  remits,  however,  were  contradicted  by  the  eiqperiments  of  Lakeman 
and  Groosmullcr  (1928). 

After  another  decade  the  subject  was  again  reformulated,  van  Cittert  (1938) 
showed  that  in  a  plane  illuminated  by  an  incoherent,  nearly  monochromatic 
source  the  optical  disturbance  is  normally  distributed.  Three  theorems  con¬ 
cerning  his  "komplex  korrelation"  were  also  given  in  this  paper  The  best 
known  of  these  theorems  expresses  the  correlation  in  the  Illuminated  plane  in 
terms  of  the  intensity  distribution  across  the  source  plane  The  analysis  of 
thif  oaper  is  in  terms  of  ensemble  averages,  but  by  invoking  the  ergodic  hypothesis, 
it  can  be  sho’vn  that  when  the  treatment  is  applicable  (Incoherent,  quasi -mono- 
ci.romatic  sources)  the  roi  relation  function  introduced  by  van  Cittert  is  equal 
to  the  zero  ordinate  of  the  complex  degree  of  coherence, 

y,  ^  yi2(0)  (1.1.2) 

A  significant  augmentation  of  tiic  theory  of  partial  coherence  was  given  bv 
'"'ernike  (1930).  In  ibis  paper  tlic  treatment  of  the  subject  is  tightly  bound  to 
t*'e  interpretation  of  Young's  interference  experiment  in  terms  of  Michelson's 
visibility.  In  fact  the  uegree  of  coherence  between  the  disturbances  at  two 


is  defined  in  that  paper  as  the  vlsibiMtv  of  ho  frirjci  jr-teined  by  afi owing  the 
light  from  these  points  to  interfere  in  i  8ui»^ble  e.  feriaient  (i.  e.  ffhc»rt  path 
differences  and  equal  intensities) .  So  valuable  has  ills  experitnent&i  defimUon 
proved  in  understanding  the  physical  aspects  ot  partial  coherence,  tnal  many 
later  authors  overlook  the  fact  that  Zernike  formulated  the  subject  analytically. 

Zeruike's  fornulation  is  applicable  to  quasi-monochromatic  fields  pro¬ 
duced  bj’  any  sov  ct  (coherent,  partially  c:^erent,  or  incoherent)  and  in  this  oense 
is  the  most  general  tff  aiment  prior  to  the  work  of  Blanc-Lapierre  and  Dumcntet 
(1955),  and  W;  l?  f^<*55)  to  be  tiscusaed  below.  The  fundamental  quantity  in  his 
analyeic*  is  ti^s  nm n  irucnaiiy  fun  Uon,  defined  as  the  time  averaged  product 
of  the  disturbcncp  ;  t  one  jKdnt  with  the  complex  conjugate  of  the  disturbance  at 
the  second  Mr.t.  The  de(?rf  ?  .if  .  oherence  is  the  mutual  Intensity  function 
suitably  nt^ui^uiaed.  Zernike's  degree  of  coherence,  v  ,  may  be  shown  to  be 

ft 

equal  to  the  zero  ordinate  (u  the  complex  degree  of  coherence  used  here, 

=  ri2(0)  •  (i-i-S) 

In  the  same  pai^er  an  a|:proxlmate  law  for  the  propagation  of  the  mutual 
Intensity  wwj  also  preseiked;  and  as  a  consequence  of  this  law  a  theorem  re¬ 
lating  to  the  mutual  Intettsity  on  a  plane  illuminated  by  an  Incoherent  plane 
source  was  determined.  This  theorem  Is  by  virtue  of  the  crgodic  hypothesis 
the  same  as  the  theorem  of  van  Clittert's  mentioned  earlier  and  is  now  terniod 
the  van  Clttert-Zemike  theorem  (of.  Bom,  M.  and  Wolf,  E.  1059,  p.  507) . 

In  1951  H.  H  Hopkins  reformulated  the  theory  of  parUal  coherence.  In 
this  treatment  the  complex  degree  of  coherence  is  defLned  In  tcrnie  of  an  Integral 
over  the  primary  sourc*..  of  the  radiation,  assumed  always  to  be  incoherent.  While 
the  arguments  of  this  paper  have  been  the  subject  of  co.naiderable  oiscusslon. 
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(of.  Wolf  (1954,  1958),  Hopkins  (1956),  Zucker  (1937)),  the  tecfmiquea  introduced 
liave  proved  very  powerful  for  many  practical  optical  problems.  When  the  formu¬ 
lation  ol  this  paper  is  applicable  (i.e.  incoherent  quasi-monochromatic  sources) 
it  is  equivalent  to  that  of  Zernikc.  Hence 

’  (l.i.4) 

wnere  is  the  degree  of  coherence  as  defined  "oy  H.  H.  Hopkins. 

In  spite  oi  tlie  uweiuiness  of  tho  results  presented  in  tiie  papers  described 
above  we  mast  point  out  three  unsatlafactory  aspects  of  their  formulations  : 

1.)  the  are  applicable  only  to  quasi-monochromatic  fields;  2.)  apart 

from  Zermke’d  wvirk,  the  analyses  are  applicable  orJy  when  the  source  is  inco- 
herer.t;  3.)  apart  frojr.  von  Laue,  each  of  the  auove  formulaUous  are  in  terms 
of  complex  functions  the  significance  of  which  is  obscure. 

The  first  of  these  cone ide rations  is  perhaps  Uie  most  significant;  for  the 
restriction  to  incoherent  sources  is  removed  by  Zernlke's  fornmlation;  and  the 
ambiguity  as  to  complex  representation  may  be  removed  either  by  dealing  ex- 
ciusivelj-’  with  leal  functions  or  by  carefully  defining  the  complex  functions.  The 
resixlction  to  quasi-monochromatic  light,  however,  is  not  simply  removable,  !n 
fact  *Jie  theory  of  partial  coherence  could  be  exte.nded  to  fields  of  arbitrary 
spectral  width  only  after  tho  introduction  of  the  cross-correlation  function  and 
»  mfore  detailed  conelderation  of  the  statistical  aspects  of  the  subject. 

These  RhOi  tcoinings  were  eliminated  in  the  formulation  of  the  subject  by  Wolf  (1955) 
iUidtliatof  Blanc- Lapierre  and  Duraontet  (1955),  Doth  of  these  formulations  are  rig- 
orouily  applicable  to  poiychri>matlc  fields  created  by  any  type  of  source  (coherent, 
partially  coherent,  or  L’scohei'ent) ;  and  both  define  tire  degree  of  coherence  in  terms  cf 
the  crt.'se-corrsiation  of  the  disturbance  at  two  {oUits  la  tho  field.  Tire  oeoentiai  JlffereiKe 
between  ‘iiees  icrmclatioiis  is  the  fact  that  W'oll  treats  the  subject  in  terms  of  carefully 


B 


! 


! 


defined  complex  functions,  while  Blanc -Lapierre  and  Dumontet  deal  in  tho  main 
with  real  functions. 

These  two  treatments  are  rigorous  and  general,  but  the  several  advantages 
of  the  complex  representation  (discussed  by  Born  and  Wolf  (1959)  and  Parrent 
(1959)  make  it  more  suitable  for  an  analysis  in  which  the  usual  optical  theorems 
for  natural  light  are  to  be  regarded  as  limiting  forms,  'Jfhe  degree  of  coherence 
as  defined  in  the  real  function  treatment,  ^  £j(t)  is  simply  the  real  part  of 
the  complex  degree  of  coherence,  i.o.. 


'BT.  4- 


ft.  1.5) 


r'or  the  reasons  stated  above  the  formulatlcn  of  coheience  tlieory  due  to 
Wolf  will  bf  ujed  exclusively  in  this  thesis.  In  the  next  section  this  formulation 
is  described  in  detail. 

Rsview  of  the  Gensrai  Formulation  of  Coherence  Theory 
As  originally  Introduces  by  Wolf  (1054)  the  iheory  of  partial  cohere  .nee  is 
formulated  in  terms  of  the  electromagnetic  fiell  The  basic  eniities  In  this  for¬ 
mulation  are  corralatlDn  matrices,  the  elements  of  which  arc  U>e  cross-curreia- 
tlons  of  the  Cartesian  components  of  the  electromagnetic  field  vectors  in  an 
appi'opriate  complex  :epre8entation.  While  such  a  geiier.'v  treatment  may  be 
necesiary  for  th;?  descriptlca  of  partai  polarlctadon  phenomena,  a  great  many 
optical  phetjornena  are  adequately  deb».ribed  by  a  scalar  wave  function. 

There  are  av  least  two  ajwroach  a  to  the  justification  of  the  use  of  a  .v’alar 
tht'ot7,  ajiart  fiom  the  exper'mei.  al  fact  that  it  Is  in  (excellent  agreement  with 
experience.  The  customary  approach  Is  to  tal.r  as  the  basic  physical  quantity  a 
allude  Cartesian  component  of  the  elv-ctric  vector,  reservirg  the  right  to  Include 
the  additional  components  if  necessary  to  dcscrib*’  some  particular  phenomenon 
(e.g.  polarization  effects). 
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An  alternative  approach  is  that  adopted  by  Green  and  Wolf  (1953).  In  this 
and  subsequent  papers  by  Wolf  (1959)  and  Roman  (1955,  1959)  a  rather  fruitful 
attempt  to  completely  describe  optical  phenomena  in  terms  of  a  complex  scalar 
function  has  been  made.  This  theory  is  not  as  yet  complete  and  hence  will  not 
be  used  here.  In  this  thesis  we  adopt  the  first  approach  and  represent  by  V^(t) 
a  single  (real)  Cartesian  component  of  the  electric  vector. 

Since  the  concept  of  Gabor's  (1946)  analytic  signal  is  fundamental  *^^0  the 
understanding  of  the  definitions  to  be  given  here,  it  will  prove  useful  tc  review 
briefly  the  meUiod  of  associating  an  aaolytic  .signal  with  a  given  real  function 
in  spite  of  the  fact  that  Chapter  2  is  devoted  largely  to  analytic  signals  and 
their  a.soociatGd  Hilbert  transforms. 

As  will  become  clear  in  the  following  development,  the  analytic  signal  may 
be  obtained  by  a  simple  generalization  of  the  method  of  associating  a  complex 
exponential  with  a  simple  periodic  ''unction.  Let  V^(t)  be  a  ’■eal  function  pcs- 
sesslng  a  Fourier  representation, 

-  f  a(c)  cos  j  ip(i')  -  27Tyt  I  di-  .  (1.2.1) 

1 

We  associate  with  V  (t)  another  real  function,  V  (t),  obtained  by  changing  the 
phase  of  each  sptict^  a)  component  of  (t)  by  n/2.  ‘fhiis 

(t)  -  f  a(c)  sin  [  -  27rp't  ]  dt'  .  (1.2.2) 

o 

The  analytic  signal,  V(t),  may  then  be  defined  as 

V(t)  -  V^(t)  +  iV^t)  .  (1.2.3) 

The  fundamental  quantity  in  this  study  is  the  mutual  coherence  function, 

I-  i  xi 

I  12'^  '  which  is  defined  is  the  complex  CiOes-correlation  of  the  analytic  | 

I 
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signals  Vj(t)  and  VjW  associated  with  the  real  disturbance  Vj* **^(t)  and  V2'^(t)  at 
two  typical  points,  Pj  and  Pg,  in  the  field,  1.  e., 

PjaW  =  (Vj(t  +  t)  V2*(t)  )  ,  (1.2.4) 

where  the  sharp  brackets  denote  time  average.*  Three  functions  derivable 
from  sufficient  importance  to  have  received  separate  names. 

The  self  coherence  function,  ^  1  ^  ^  conaple*  auto-correlation 

of  the  analytic  signal  asooclated  with  the  disturbance  at  a  typical  point,  P^,  in 
the  optical  field, 

rn(T)  »  Ij(r)  =  (Vj(t  +  t)  Vj*(t)  ).  (1.2.5) 

As  will  be  proved  in  a  later  chapter,  the  intensity  at  a  typical  point  in  the 
field  is  the  zero  ordinate  of  the  self  coherence  function,  i.  e,, 

rii(o)  s  ij(o)  .  (yj(t)  v/(t))  .  (1.2.6) 

This  result  may  be  regarded  as  a  generalization  of  the  familiar  theorem  for 
monochromatic  light. 

The  Fourier  transforms  of  the  mutual  coherence  function  and  the  self 

coherence  function  are  termed  the  mutual  spectral  density, 

^  ♦  ♦ 

spectral  density,  ^  respectively;  l.e,, 

oo 

'  (1.2.7) 

:=  0  <  0  , 

*  The  precise  form  of  ‘hie  tlm-®  average  is  discussed  In  Chai^er  2. 

**  It  will  be  proved  in  Chapter  2  that  the  cross-correla(lon  of  two  analytic 
signals  la  itself  an  analytic  signal,  and  hence  Raf*'/  contains  only 
positive  frequency. 


and 
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fjlM  -  J r„(r)  6=*”"'^  dr 

»oe 
-  0 


I'  >  0 

<  0  , 


(1.2.8) 


In  terms  of  these  functions  the  complex  degree  of  coherence  (function)  is 
defined  as  the  normalized  mutual  coherence  function,  Miere  the  normalization 
factor  is  the  square  root  of  the  product  of  the  intensities  at  and  P2; 


vTii^nv^ 


(1.2.9) 


By  appealing  to  the  Schwarz  inequality,  it  can  be  shov^’n  that  the  modulus 
of  bounded  by  zero  and  one, 

0  ^  I  yi2(T)|  :L  1  .  (1.2.10) 

and  these  extreme  values  characterize  by  definition  complete  coherence  and 
complete  Incoherence,  respectively.  It  nay  be  shown  (see  Wolf,  1955)  that 

^  identified  with  the  visibility  of  the  fringes  obtained  by  causing 
Vj*^(t)  and  V2^(t)  to  interfere  with  a  path  difference  of  vt,  where  v  is  the 
velocity  of  light  in  the  medium,  assumed  to  be  homogeneous.  It  is,  therefore, 
clear  that  the  above  definitions  are  in  agreement  with  the  conaideratlon  mat 
coherent  light  interferes  and  incoherent  light  does  noi.  In  fact  the  definitions 
used  here  may  be  regarded  as  a  rigorous  generalization  of  those  Introduced  by 
Zernlke  (1938). 

Tlie  essential  mathematical  structure  of  coherence  theory  Is  contained  in 


equations  (i.2.1)  through  (1.2.10)  except  that  Ihe  form  of  the  time  average  has 
not  been  explicitly  given.  'Oils  omission  will  be  remedied  In  Chapter  i  wiiere 
alternative  definitions  of  the  time  average  are  discussed. 
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CHAPTER  2 

MATHEMATICAL  DEVELOPMENTS 


I 

I 


The  definitions  Introduced  in  the  previous  chapter  employ  the  concept 
of  Gabor's  analytic  signals  (see  Gabor  (1C46));  and  since  the  development 
of  subsequent  chapters  makes  extensive  use  of  the  properties  of  analytic 
signals,  it  will  prove  helpful  to  discuss  in  some  detail  their  mathematical 
properties  at  this  point. 

The  method  given  in  Chapter  1  for  assjciating  an  analytic  signal  with 
a  given  real  function  is  patterned  after  the  introduction  by  Wolf  (1955)  and 
is  useful  in  understandir^  the  significance  of  the  analytic  signals.  However, 
for  the  purpose  of  the  analysis  given  here  it  is  more  convenient  to  IntroAice 
an  alternative  definition  in  terms  of  Hilbert  transforms.  The  equivalence  of 
these  two  definitions  will  be  demonstrated. 

In  Section  2.2  several  convolution-type  theorems  for  analytic  signals 
are  given.  Section  2.3  presents  a  numbt^r  of  theorems  concerning  the  cross- 
correlation  of  analytic  signals.  In  the  final  section  the  form  of  the  most 
g^^eral  unimodular  anai^/tic  signal  is  given. 

2.1  Hilbert  Transforms 

The  Hilbert  transform,  v\t)  .  of  the  real  function  (t)  ij  defined 
by  the  relation 


V^t) 


and  V^(t) 


Here  the  Integral  is  the  Cauchy  principal  value  integral  which  may  be  de 
fined  as 
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. (2.1.2) 

However,  in  almost  every  instance  it  wih  prove  more  convenient  to  evaluate 
the  principal  value  integral  by  contour  integration.  To  rewrite  (2.1.2)  in 
terms  of  contour  integrals  we  replace  t  by  a,  (z  =  x  +  iy  where  x  and  y 

are  real),  and  assume  that  the  function  V**  (z')  is  continuous  in  some  neigh¬ 
bourhood  of  e'  =  E. 

The  principal  value  Integral  may  hen  be  defined  by 


where  th '  integrals  on  the  right  are  line  integrals  along  the  open  curves 
and  C  .  The  curves  and  C_  extend  along  the  real  axis  from  left  to 
right  curvtr^  in  semi-circular  arcs  respectively  above  and  below  the  point 

i'  =  T  as  in  Figure  1.  The  value  of  the  line  Integrals  in  (2.1.3)  are  of  course 
independent  of  the  actual  path  of  the  curved  part  of  C ^  and  C_  as  long  as  they 
do  not  Include  singularities.  The  line  integrals,  therefore,  r«mab>  unchanged 
in  the  limit  of  vanishing  radius  for  the  seral-clrcular  arcs.  The  continuity 
of  V**  (e')  in  the  neighbourhood  of  e*  =  e  assures  that  In  the  limit  the 
contributions  from  the  \xpper  and  lower  arcs  will  cancel  each  other;  thus  the 
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C+  To  lofinit/ 


C-  To  In-finilv 


I  iy  FIGUftff  1. 

equivalence  of  the  definitions  of  principle  value  is  den'.onairated. 

The  integrals  in  (2  1.3)  are  line  Integrals;  the  patho  by  which  and  C 
are  closed  will  of  course  depend  on  the  form  of  {z')  in  any  given  problem 
Ln  terma  of  the  Hilbert  tran8fo.’*m  the  analytic  signal  may  be  defined  as 
follows:  letV^(t)  be  a  real  function  of  t  such  that  its  Hilbert  transform, 


vVt),  exists.** 


The  analytic  signal  may  then  he  defiricd  as 


V(t)  =  v‘’(t)  +  ivV) 


(2.1.^) 


Actuallv  (2.1.3)  is  a  more  general  definition  Includiiig  (2.1.2)  as  a  special 
case;  for  (2.1.'^;  will  not  exist  if  v’^(z’)  has  a  pole  at  z'  =  z  while  (2.1.3) 
will  always  exist  except  foi  tiie  cases  of  an  essential  singularity  or  btanch 
!X)int  at  z'  =  z. 

**  A  sufficient  condition  for  the  existence  of  tl.,  Hilbert  transform  is  tl>at 
V^z)  be  square  Integrable.  Howevc:,  as  with  the  Fourier  transform  the 
necessary  conditions  for  the  existence  of  the  Hilbert  transform  are  nc. 
known.  Since  many  of  the  functions  oc'^urrlng  in  thlo  chapter  are  not  square 
Integrahle  but  do  posess  Hilbert  transforms,  we  assume  the  weaker  condi¬ 
tions  that  the  transform  exists  and  that  the  inversion  theorem  is  valid 


\ 


>*l 
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The  equivalence  of  this  and  the  previous  definitiOHj  (1.2. i)  -  (1  ?.3), 
niay  be  demon&traied  as  foliovvs:  express  V*  (t)  as  die  Fourier  integral 


i' t )  =  j  3  (v)  cos  [  ^  ( i-')  -  2  ji  ]  d  V 

0 


(2.1.5) 


Take  the  Hilbert  transform  of  Doth  sides  of  (2.1.5)  after  setting  1  =  z  and 
exi:re83nig  ihe  cosine  in  terms  of  complex  exponentials.  Thus 


00  00 


v^(?)  -  —  /  a(c)  f  2 - - 

0.1T  f  c'  - 


J~  J  [0  ('■')  -  ] 


dz  ’di/ 


<3C 

+  —  /  a(i')  i  - 

2jr  J  J 


i  [  ^  (i')  •'  2  ffi'z'  ] 


(2.1.6) 


J  z'  ~  z 


dz’dc 


♦  H 


Ihs  order  of  integration  ha>«  Deen  interchanged  and  the  inner  integrals  are 
intorpr..tsft  cs  die  line  iclegraU  in  Figure  i,  'D»e  contour  is  closed  at  in¬ 
finity  below  the  axis  in  the  first  Inlsgral  and  at  infinltj'  above  the  axis  in  the 
second  integf'al.  Since  there  are  no  poje.s  within  the  contour  in  either  inte¬ 
grand  except  at  z'  =  z,  we  obtain  by  the  residue  theorem 


V  (t) 


a(v)  sin  r  0  ( i^)  2iTct  ]  dc  ,  (2.1.7) 


r*-  >f 

>N 


Wiiere  we  have  replaced  z  by  t  after  performing  the  integration.  Comparison 
of  (2.1.7)  and  (1.2.2)  establishes  the  equivalence  of  the  two  definitions. 

It  will  be  convenient  for  later  deveiopra  ?nt8  to  list  here  several  theorems 
concerning  Hilbert  trarsforms  (c.f.  Bateman  Monuscrlpt  Project  1934). 
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Table  1 

Theorem 

I 

n 

ra 

IV 

y 

2, 2  Ccnvolutlors  ItliieoreniS  for  Arm  lytic  Slgtmla 

We  1  rttiulre  feveral  theorems  cojicernlrig  the  cross-correifttlon  of 
ar-alytic  Since  these  ^heorsn*.*  cio  not  seem  to  appear  Ui  the  literature, 

we  ehall  gi^o  a  detailed  derivation  of  them  ai  this  potr-t  Becfeu?s  of  tht^  com¬ 
plication  LitroJtuced  by  the  lime  average  in  the  cross-correlatlcn  hinctlOii,  we 
shall  first  viemoriStrate  these  theorejns  in  *erm«  of  convolution  integrals  anu 
i  lat^r  derive  the  ccrrelatiun  thooxems  from  them. 

'  Theorem  VI :  'Hw  convolution  ui  ‘uro  reel  fimctions  fj(t)  arid  {^(t)  ie 

j  equal  ilie  cOiSVCiytSon  of  their  Hilbert  tran-sfoisiis,  g|(t)  sod  gjtt}  0^  the  asme 
;  order),  l.e., 

t 

f 

1 

1 


f(t) 

1  /  f(t’)dt' 

f  J  t'  -  t 
<00 

V^(t) 

y\t) 

-  V^(t) 

V^(t  +  a) 

V^(t  +  ; 

V^(t  at) 

tvV+  at) 

+  a)  V^t) 

(t  +  a)  V^(t) 

-i  Vd) 

dt 

d  v^(t) 
dt 

V*’(t)dt 


K 


r 


00 


-«Q 
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dt 


+  T)  dt  (2.2.1) 


-« 


The  proof  of  this  theorem  can  be  obtained  with  the  help  of  Theorems  I  and  n 
of  Tabic  1.  Let 


?(t)  [  fj(t)  fgCt  +  T)  dt  . 


(2.2.2) 


Appljing  the  invereion  theorem  (Theorem  I,  Table  1)  twice  to  both  sides  of 
(2  2.2),  we  obtain 


IT 


fj(t)  f2(t  ^  T)Qt  dTdT'  . 


.  .  ^2.2.3) 


Usliv  'Hieorem  U  of  Table  1  after  inverting  the  order  of  integration,  we  find 
that 


^ — -  Ai(0  *2'^ 

Jte  *Ki 


(2.2.4) 


Jn^iod.iCiM  \  nev'  variable  t  -  t  + 


W9  obtain 


F( '") 


^^(t) 

t  -  (t  -  T-) 


(2.2.5) 


\ 


18 


Uslns  Xlieoi^m  E  again  we  ^innliy  hp  ve  the  rseidt  that 

•C 

-  J  -  T")  ^ifidr .  (2.3r«) 

On  setting  T  -  r'  =  t  (2.2.8)  reduces  to  (3.2.  j)  aiid  the  shfHorem  is 
demonstrated. 

Two  special  cases  of  theorem  (2.2, i)  are  o!  nartlcuiar  iotereat.  The 
first  case,  £^(1)  -  V^(t}  (a  i,2),  leads  to 

*c  aa 

J  Vj'’(t)  V\J{t 

^  'tl  -4515 


T)dl  = 


v,'(t)  Vj^Ct  f  T)dt,  (2.2.7) 


Toe  second  case,  L  (?)  =  v^^(t)  an^.  f2(t)  =  icaos  to 


j  +  T)dt  r  -  J  Vj’(t)  Vg^'Cl  +  Tidt.  (2.2.8) 


r 


Using  (3,3  7)  and  (3.2.8),  we  ma>  i>b'aia  the  convolution  thsejrtms  foi  aaalytlc 
c  tgnak  from  which  some  o,1  the  advantages  of  this  comivia*  ?ormu?stiojj  will 
be  ciesiv  Consider  the  inlegml 


f  Vjd)  7,‘(t  ^  r].ii  . 

J  ^  - 


[Xt.j] 


wy  .'reV^it)  (#  ^  1,2)  is  of  the  form  of  (2. 1.4).  SybaUtutix^  irom  (2. 1.4) 
into  (2.2.9)  and  using  (2.2.7)  and  (2.2.3/  obtain 


t 

ii 
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«0  CO 

J  V^a.)  it  +  T)dt  =  2  j  Vj^(t)  V2*‘(t  +  T)dt  (2.2.10) 

-00  -CO 

00 

+  2i  J  Vj^(t)  V2^(t  +  T)dt  . 

-00 


Equation  (2.2.10)  <?xp£«sses  tha  usefol  result  that  the  real  pane  of  the  convolu¬ 
tion  of  two  analytic  signals  is,  apart  from  a  factor  of  two,  tlie  convolution  i.f 
the  real  fur  ‘ont;  w  l^h  which  they  are  associated.  Another  useful  property  of 
ih<5  anal  t- si^msls  is  seen  by  putting  Vj(t)  equal  to  V2(t)  in  (2.2.10).  Thus 


? 

i  I  Vj(t:.  Vj*(t  *  T)dt  -  j  Vj^(t)  v/(t  +  7)dt 


•<00 


(2,2. il) 


+  i  J  Vj^'d)  v/  (t  *  r)dt. 

-oo 

We  may  now  evaluate  the  second  integral  on  the  right  as  follows;  using  (2.2.8), 

CO  QO 

J v/(t)  VjVt  +  r)dt  =  -  J V, V)  Vis’ll  +  T)dt  , 


or  evaluating  (2, 2. 13)  at  t  -- 0  we  find 
00 

J  Vj'(t}  v/d)(lc  ^  0  , 


.  (2.2.12) 


(2.2.13) 


—  1 

l.e.,  Vj  ti)  and  Vj  (t)  are  orthogonal.  Evaluating  (2.2. 11)  at  t  =  0  and 
using  (2.2.13),  we  obtain 


(2.2.14) 


/  Vj(t)  V,  (t)d(  =  ?.  /  V/(t)  Vj^(i)dt  . 


i, 

I 

H 


* 


4- 
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Equation  (2  2,14)  expi^f  ses  a  second  important  property  of  the  analytic  signals; 
namely,  the  integi'ai  over- all  time  of  the  squared  modulus  of  an  analytic  signal 
is,  apart  from  a  factor  of  two,  the  Integral  over- all  time  of  the  square  of  the 
real  function  with  whic't  it  ia  associated. 

An  additional  theorem  on  the  convolution  of  analytic  8i;;,nal6,  which  will 
pro\i;  important  in  the  next  section  is  : 

Theorem  Vn  :  The  convolution  of  two  analytic  signals  is  itself  an 
analytic  signal.  To  demonstrate  this  theorem  we  write 


F(r)  -  2  V2*(t  -  T)dt  . 

Using  (2.2.7)  and  (2.2,8),  (2.2.15)  can  be  rewritten  ae 

OC  00 

Fit)  =  J  Vj^'{t)  V./H  4-T)dt  +  yv,^(t)  +  7)dt 


(2.2.15) 


-•c 


i  i 


/  *  * 

I  VjV  +  r)dt  -  r  Vj‘(t)  V3*’(t+  T)dt|. 


.  (2.2.16) 

Denoting  by  R  and  I  the  real  and  imaginary  parts  respectively  and  taking 
the  Hilbert  *"ansform  of  the  real  part  of  F(t)  we  obtain 


«0 

1  I 

^  j  r  -  T 


df 


-  I  {f  (t)], 


(2.2.17) 


and  'niecrem  Vn  is  demonstrated. 

The  results  established  Ln  thic  section  are,  of  course,  applicable  only 
to  convolutions  and  the  definitions  Introduced  in  Chapter  1  are  In  terms  of  correla- 
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tlons.  We  shall  shOAw  in  the  next  sertion  that  these  theorems  are  also  valid 
for  the  cross-correlation  of  analytic  signals  if  tlie  time  average  is  suitably 
deflred. 

2.3  The  Cross-Correlation  of  Analytic  Slgnale 

The  complex  cross-correlation  fimction  may  be  defined  in  several  ways. 
The  cusiomary  definition  ia  (c.f.  Daverport  and  Root  (1958),  page  70) 

T 

^j2(t)  -  lim  f  fj(t)  f2*(t  ^  r)dt.  (2.3.1) 

-4 

However,  for  our  present  purpose  the  mutual  coherence  fiaiction^  ’ 

defined  as 

SO 

ri2{^  -  lira  3Y  rVi(T,t  +  T)  V2*(T,t)dt,  (2.3.2) 


where 


y(T.t)  = 

v'‘(T,t)  - 

l.pV‘ft) 

(2.3.3) 

ami 

(  V*'(t) 

hi 

<  T  ^ 

=  ( 

1  1 

} 

(2.3.4) 

1  ^ 

hi 

>  T  j 

and 

‘9- 

TV‘(t) 

-7 

V*'(T.t')dt' 
t^  -  t 

(2.3.5) 

The 

different  notation, 

v''(T,t)  and 

jVh),  is 

used  smee, 

while  /(T,t) 

vanishes  for  1 1  |  >  T,  Its  Hilbert  transform  ,j.V^  (t)  will  no*  in  general 

H 
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1;'^  ,  . 


vanish  in  this  range.  Since  V^(t)  is  assumed  to  be  everywhere  finite,  the 
function  V*^(T,t)  is  square  integrable:  its  Hilbert  transform  ^V^(t),  is 
therefore  also  square  integrable  (c.f.  Titchmarsh  (1948) ).  Consequently, 
all  the  required  Hilbert  transforms  exist.  Tlie  equivalence  of  the  two  defini¬ 
tions  of  an  average,  of  the  type  (2.3.1)  and  (2.3.2)  is  discussed  by  Born  and 
V/olf  (see  Born  and  Wolf  (i960) ,  page  496).  For  our  purposes  it  is  more 
convenient  to  employ  the  definition  (2.3.2)  and  accordingly  the  sharp  brackets 
are  defined  as 


(v/(t+  T)  Vj  (t) 


)  =  hiL  / 

■P**  J 


Vj(T,t  +t)  Vj  (T,t)dt. 


The  presence  of  the  parameter  T,  In  V(T,t)  in  no  way  affects  the  argu¬ 
ments  of  the  previous  section;  and  the  operation  will  commute  with 

the  integrations  involved,  since  T  and  r  are  independent  variables.  The  theorems 
established  in  Section  2.2  may,  therefore,  be  bdeen  over  mutatls  mutandis  for 
the  cros  -correlation  function  (^2(t)  .  Thus  If  we  adopt  the  notation 


<  v/(l  .  7)  V/(t)  ) 


We  mzy  summarize  the  principal  theorems  as  follows  : 

ri2(^  =  2  [  rj2‘''’(T)  +  1  r 


(2.3.6) 


(2.0.V) 


r,2"(T)  7  1  r;/‘(r)  ' 

. 

(2.3.8) 

0  • 

(2.3.9) 
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Setting  V|(t)  equal  to  V 2 (t)  and  r  =  0  in  (2.3.8)  and  using  (2.3.9),  we 
obtain  the  result  anticipated  in  Chapter  1,  Section  2,  viz., 


rii(0)  =  2  rii'*^(0)  =  2  I(Pj) 


(2.3.10) 


where  I(Pj)  is  the  intensity  at  Pj. 
the  arguments  of  this  section  tnat 

rj2^^T)  =  -  r'j2^’^(r) 


Further  It  follows  from  (2.2.19)  and 


1  f  4/^r')dT' 

J  r'  -  1 

-00 


(2.3.11) 


Thus  the  mutual  coherence  function  is  an  analytic  signal. 

For  the  sake  of  continuity  in  latei  arguments  we  include  at  this  point 
two  lemmas  concerning  complex  cross-correlation  functions.  Both  of 
these  lemmas  .re  well-known  and  follow  immediately  from  the  definition 
of  and  the  atatlonai  lty  condition ;  they  are,  therefore,  given  here 

without  proof: 

Lemma  I 

r2i(-T)  =  ri2V)  .  (2.3.12) 

and 

Lemma  II  ,, 

■  (2.3.13) 

The  hooked  notation  is  used  to  denote  temporal  Fourier  transforms,  i.e., 


f 


J 


2xrlr  r 


dr 


>  0 


=  0 


r  <  0 


(2.3.14) 


i 


[ 
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The  spectrum  of  r*j2(‘^  1®  zero  for  half  the  frequency  range,  since  the 

mutual  cohei*ence  function  Is  an  analytic  signal. 

2.4  The  Iitost  General  Unimodular  Analytic  Shiml 

We  shall  be  concerned  in  Chapter  3  with  the  determination  of  i  mitlng 
forms  of  the  mutual  coherence  function  fcr  coherent  and  incoherent  fields. 
The  form  of  172^^)  for  coherent  fields  will  be  shewn  to  be 


(2.4.1) 


1^12^ 


where  e  is  a  unimodular  analytic  signal.  The  determination  of 

involves  the  solution  of  a  singular  integral  equation;  and  since  this 
solution  is  rather  lengthy  and  purely  mathematical  in  nature  it  will  be  included 


in  this  chapter 


‘♦12(^1 


Since  e  is  an  analytic  signal,  its  real  and  imaginary  parts  are 

Hilbert  transforms,  i.e., 

«o 

)  dr' 


sin  ^i2('r) 


1  / 

^  J  t'  -  T 


(2.4.2) 


Using  the  inversion  theorem,  Theorem  I,  Table  1,  we  njay  write 

T')dT' 


cos  = 

Combining  (2.4.2)  and  (2.4.3)  we  obtain 


i  r ^12^^ 

^  J  -  r 


(2.4 


A 

’  J  r 


dr 


(2.4.4) 
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Here  we  have  dropped  the  subscript  l,  2  since  the  development  of  this  section 
is  not  concerned  with  the  space  dependence  of  ' 

Sqvstlon  (2.4.4)  is  an  integral  equation  for  and  may  be  recog¬ 

nized  as  a  singular  form  of  a  Carleman-type  equation,  i.e., 


f(T)  = 


OQ 

+  X(T)  , 


(2.4.5) 


where  x  is  a  generally  complex  constant.  The  general  solution  to  (2.1.5) 


is  (c.f.  Tricomi  (1957),  page  175) 


'  rrV?  ^  I 


(2.46) 


However,  since  in  (2.4.4)  X  =  l/v  and  x(t)  =  0,  the  solution,  (2.4.6) 
becomes  in  this  case  indeterminate.  In  fact  it  is  clear  that  (2.4.5)  can  have 
no  S4.  ution  for  a(T)  =  0,  except  possibly  for  the  eigenvalues,  X  =  +  i/v  . 
This  singular  form  is  solved  here  by  utilizing  the  pro]>ertle8  of  analytic  func¬ 
tions.  To  this  end  wc  replace  the  real  variable  r  by  a  complex  variable, 
z  X  +  iy;  and  write  (2.4.4)  as 


f(*)  ^  ^ 


f(z’) 

z  ’  -  z 


(2.4.7) 


where 


f(z)  . 

Since  ^  (x)  is  a  real  function,  we  have 


(2.4.8) 


^(i)  ; 


(2.4.9) 
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here  and  throughout  tlie  rest  of  this  section  we  use  the  bar  to  denote  complex 
conjugate.  In  tiie  subsequent  discussion  we  shall  treat  only  the  case  X  -  -i/v 
since  the  argument  is  easentirlly  the  same  for  the  positive  eigenvalue. 

It  has  been  pointed  out  in  Chapter  1  that  the  term  analytic  signals  derives 
from  the  fact  that  these  functions  are,  when  considered  as  a  fwction  of  a  com¬ 
plex  variable,  analytic  in  half  the  complex  plane.  Since  extensive  use  of  this 
property  is  made  in  this  section  we  shall  digress  briefly  and  demonstrste  it. 

Let  U(t)  be  any  function  of  the  real  variable  t  such  that  the  integral 


F(z) 


(2.4.  iO) 


exists.  Then  the  function  F(£)  is  analytic  in  the  lower  half  of  the  complex 
plane  (c.f.  Whittaker  and  Watson  (lOfK}}  page  92).  Equation  (2.4.10)  may  be 
rewritten  as 


Hi  U(2) 


where  Cauc.sy's  theorem  has  been  used. 

In  the  limit  as  z-*x  from  below  the  real  axis 


(2.4.11) 


«o 


=  u(x)  -  1  j  .  (2.4.12) 


From  (2.4.12)  It  Is  clear  that  the  analytic  signals  may  be  regarded  as  the 
limit  as  the  real  axis  is  approached  of  a  function  analytic  in  half  the  complex 
plane.  Oi,  conversely,  if  in  an  analytic  signal  the  real  variable  is  replaced 


I 
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by  a  complax  one  the  resulting  function  Is  analytic  in  haU  the  complex  plane. 

For  our  present  problem,  the  jolution  of  (2.4.7),  the  domain  of  analyticity 
of  f(e)  may  be  extended  to  include  tne  real  axis  since  1.)  (2.4.9)  excludes 
the  fossiblllty  of  poles  or  essential  singularities  on  the  real  axil,  and  2.) 
physically  we  require  f(t)  to  be  unique  which  excludes  branch  points  on  the  real 
axis. 

Apart  from  the  trivial  solution  of  f(z)  identically  eqtul  to  zero,  the 
function  f(x)  can  have  at  most  Isolated  zeros  in  the  lower  half  plane.  This 
conclitsion  follows  from  the  analyticity  of  f(z)  in  that  domain  (c.f.  Titchmarah 
(1939)  page  88). 

Equation  (2.4.9)  expresses  the  value  of  f(z)  at  all  points  in  the  lower  half 
plane  in  terms  ot  its  values  at  conjugate  points  in  the  uf^r  half  plane.  Further, 

(2.4.9)  indicates  that  corresponding  to  every  aero  in  hte  lower  half  plane  there 
is  a  pole  at  the  conjugate  point  in  the  upper  ball  plane  and  ctmversely  zeros  in 
the  upper  half  plane  correq;)ond  to  poles  in  the  lower  half  plane.  Hence,  there 
are  no  zeros  in  the  i9per  half  plane  and  the  singularities  in  the  xipper  half  plane 
are  .solated.  Furthermore,  the  singularities  in  the  upper  half  plane  are  poles 
and  not  essential  singularities;  for  by  the  Weierstrass  theorem  (c.f.  Titchmarah 
(1939),  page  93)  in  every  neighbourhood  of  un  iaola:6d  essential  singularity  the 
function  tends  to  any  given  limit  an  infinite  number  of  times,  and  this  behsviour 
woxild  by  (2.4.9)  be  reflected  irao  the  lower  half  plane.  That  is,  an  eaci^^ntial 
singulnrity  in  the  upper  half  plane  would  necessarily  correspond  to  an  essential 
singularity  in  the  lower  half  plane.  The  preceding  argument  does  not  exclude  the 
possibility  of  essential  alngularltios  at  infinity  since  infliUty  is  not  an  isolated 
point.  In  the  same  way  a  branch  point  in  the  v^per  half  plane  Is  excluded  since 
it  would  Imply  a  branch  point  in  the  lower  half  plane  contrsdlcting  the  requirement 

► 

i 


28 

of  analyticlty.  We  m»y  conclude,  therefore,  that  f  (z)  ie  a  meromorphic  func- 
tion,  i.  e.,  ita  only  singularltied  for  finite  z  are  poles. 

By  a  modification  of  Hadamard's  theorem  (c.  f.  Tltchniarsh  (1839),  page 
284g  where  a  full  discussion  of  the  concepts  employed  below  is  given)  any 
meromorphic  fimction  f  (^)  may  be  e:iqpressed  in  the  form 


f(z)  = 


ri'  s| 

n=l 


where  Q(z)  is  a  polynomial  oi  order  N 


;  HEfe-,, 

ii*l  Pn 


the  canonical  product  of  the  primary  factors. 


E(u,p)  *  (1  -  u)  6 


u 


P 


is 


(2.4.15) 


(2.4.18) 


a^^  and  b^.  are  the  zeros  and  poles  respectively  of  f(*). 


The  genua  of  the  canonical  product  satisfiea  the  inequality 


<  p  0  =  1,2)  , 


(2.4.17) 


where  p  la  the  order  of  the  mciomorphlc  function.  The  order,  N,  of  the 
polynomial,  Q,  also  satlafiee  the  inequality,  (2.4.17),  i  e.. 


N  <  p  .  (2.4.18) 

From  the  previous  argument  It  is  clear  that  the  poles  and  aeroa  oc';ur 
at  conjugate  points,  i.  e., 


(2.4.19) 


Therefore,  using  (2  4.16),  (2.4.15)  may  be  written  as 


f(z) 


m 


"■  1  -f '  j=i 
-Ale 


ib.  +  Z)  /  1 

J  n=l  a„ 

1  1 

zl  ■ 

■  '57 

J! 

L 

"j. 

n=I \  j  ^ 


(2.4.20) 


where  we  hp.ve  set 


m 


Q(i!)  -  Z)  b.  , 

1=0  J 


(2.4.21) 


and  m  Is  the  largest  o£  U»e  integers,  pj,  Oj,  N.  After  setting 


it^  =  2:  u 


n  =  1  la 


1 


(2.4.22) 


where  Cj  Is  a  real  constant,  (2.4.20)  can  be  rewritten  as 


na 


f(z) 


/ 

{  f 

1 

’  a 

n 

>1 

1 

-JlI 

i  E  (b^  op  J 


e  j=0 


)'■ 


(2.4.23) 


However, 


Hr)  -  e 


l(b  +  c  )z 
m  m' 


m 


(2.4.24) 


when  I  *  !  — ►  *0  ;  but  (2.4.24)  has  m  poles  and  m  zeros  equally  spaced  on 
the  circle  at  infinity;  and,  therefore,  11  m  >  1,  f(a)  has  slngulariiies  In  the 
lower  half  plane.  We  conclude,  therefore,  Uiat  the  meat  gentral  allowable  solu- 
tlcn  is  with  ni  =  1.  The  expression  (2.4.22)  becomes  flnalh 


I 

I 
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a  meroEiOrpliic  function  of  order  1.  Thus  we  have  established  ^he  following; 

Theorem  Vm  :  ITje  most  general  unijnn.»ciulkf  a^iaiytic  signal  is  a 
meraxQorphic  functic?  of  order  one  with  zeros  z  only  in  the  lower  half 

plane  and  with  poles  at  conjugate  poindii:  in  the  ii|^r  half  plane  and  is  given 
by  the  formula  (2.4.25)  where  $  and  y  are  real  constants  arid  the  imaginary 
part  of  a^^  is  gieater  Uwin  zero. 
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CHAPTER  3 

SOME  IMPLICATIONS  OF  COHERENCE  AND  INCOHERENCE 

In  Chapter  2  it  was  pointed  out  that  the  modulus  of  the  degree  of  coherence 
is  bounded  by  one  and  zero  and  that  these  extremes  are  by  definition  charac¬ 
teristic  of  coherence  and  incoherence  respectively.  For  quasi-monochromatic 
light  these  definitions  are  in  accord  w.th  the  consioeration  that  coherent  light 
interfere''  and  incoherent  light  does  not.  (See  Wolf,  1954.)  However,  in  this  i 

chapter  we  arc  concerned  witli  the  implications  of  coherence  and  incoherence 
on  the  form  of  fToCT),  for  fields  of  arbitrary  spectral  width.  In  a  detailed  j 

j 

analysis  of  such  fields  in  which  the  terms  "coherent  limit"  and  "incoherent  j 

limit"  are  to  have  a  clear  and  unambiguous  meaning  the  above  defiritions  must 
be  more  precisely  stated. 

The  ambiguity  which  arises  in  the  study  of  poly ch^-oma tic  fields  s'vraa 
from  the  fact  that  the  modulus  of  the  degree  of  cnhertnce  between  the  dietur- 
bancee  at  two  points  is  a  function  not  only  of  the  position  of  th^  two  points  but 
also  of  the  time  delay  r  .  Thus  it  is  p<.:'88lble  that  for  soir.e  vilues  oi  r  , 
say  T  =  and  T  =  Tj  ,  |yi2('i)j  "  ^  ^  ° 

the  same  pair  of  disturbances  in  such  fields,  and  tne  limiting  concepts  of 
coherence  and  Incoherence  become  amblgiious. 

niis  difficulty  does  not  aiiso  in  the  study  of  quasi- monochromatic  light 
for  the  approximation  characterizing  such  fields  makes  the  modulus  of  the  de- 
greo  of  coherence  independent  of  for  all  values:  of  r  fur  wMch  the  theory  is 
applicable.  TTus  conclusion  is  evident  from  Uie  following  con3lde''ation<5.  For 
quae ‘.-monochromatic  light  the  mutual  coherence  funcllcn  ior  sufficiently  small 
'  T  j  is  given  by  (see  Appendix  2) 

V 

.  ( 


'.Ip 


I* 


=  '^2(0)  e 


-2iriPT 


Hence 


1^:/. 


-SffiX'T 


Ai' 


y,2(^  -'  y,2(o) 


«Oj-; 


,  « 


I  '  n 
'  •»  ■ 


If)’  “ .  , 


Here  P  is  the  mean  frequency  and  £u>  Le  the  spectral  wldi-^  of  the  illurnLnatlon. 

"nie  ambiguity  which  arises  in  discussing  polychroinatic  ilght  will  be 

removed  if  for  the  limiting  cases  of  coherence  ani  incoherence  we  demar>d 
that  the  modulus  of  the  degree  of  coherence  be  t- independent.  Accordingly, 
the  following  definitions  are  introduced  : 

I.  The  DiSTURBANCSH  V,  (t)  end  Vjlt)  will 
be  described  as  coherent  if  j  >'  ^(7)  i  ^  1 

!  i<3  '  ' 

for  all  r  and  inci'he rent  if  ~  ^ 

for  all  T. 

II.  An  OPTICAL  FIELD  will  be  said  to  be 

[lna>he^^ni  I  disturbances  at  all 

5>a.*rs  of  points  in  the  (teid  \ 

(incoherent' 

'.hus  i  ^  for  ?  coherent  field,  and 

|>',.,(t)|  0  for  an  Incotw! rent  field. 

Since, in  tr.e  paaLdstailed  ar.#tysl8  hae  Let-n  limited  tr  quasi-monochromatlc 
light,  it  to  clear  from  the  above  discussion  that  this  modification  of  the  defini¬ 
tions  canr».jt  lead  to  contradirtun  witii  die  wrrk  of  earlier  writers. 


The  inalysis  oi  the  general  properties  of  partially  coherent  wave  fields 
is  conveniently  performed  in  terms  of  the  mutual  coherence  function^ 
rather  than  the  oegree  of  coherence,  shall,  therefore,  examine 

the  implications  of  the  itbove  definitions  concerning  the  form  of  tiie  mutual  co¬ 
herence  function  in  coherent  and  incoherent  fields.  Tl.iese  limiting  forms  of 
r^2('^)  will  prove  useful  in  examining  tine  extremes  of  the  propagation  law 
obtained  in  Chapter  4  and  in  recovering  tiie  familiar  forms  of  elements 
optical  wave  theoiy  from  the  generalized  transfer  functions  intrrd?if,ed  iu 
Chaptei  5.  Apart  from  these  immediate  applicatk-  Tr,  the  liniitint;  foi  ms  of 
the  mutual  coherence  function  will  prove  useful  in  thci  •  own  right  since  they 
will  provide  insight  into  the  structure  and  questions  of  existence  of  coherent  and 
Cicohorent  fields. 

3.1  Wave  Equations  for  the  Propagatten  of 

The  determination  of  the  form  of  q2^^^  coherent  ani  incoherent 
fields  Will  be  seen  to  be  intimately  connected  v/ith  the  fact  tiiat  in  vacuum 
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( r)  is  pi’opaf'  (rigorously)  accord  jig  lo  the  wave  equations 


--  -  (S  1,2). 

:  0  r* 


(3.1.1) 


Here  V  is  the  Laplacian  operator  in  the  coordinates  of  P  and  c  is 

the  velocity  of  light  in  vacuum.  For  Ihis  reason  it  will  prove  helpful  to  digress 

* 

briefly  and  derive  these  equations. 

As  pointed  out  In  Chapter  1  we  consider  Uie  optical  field  to  be  characterized 
by  a  real  sc£.lar  function,  V‘(t),  which  in  vacuum  satisfies  the  equation 


These  equations  were  originally  derived  by  Wolf  (11155)  by  different 
arguments. 
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,2  vr(„  ,  i 

8  r 


(3.1.2) 


We  recall  that 


00 

=  -  f  ylJn^ 

f  ■  t 


(3.1.3) 


Operating  on  both  sides  of  (3. 1.3)  with  the  Laplacian,  V  ,  we  obtain 

00  «»'  8^  V'*(t*) 

^2  v‘(t)  =  ^  r  vV(t')dt'  ,  r  4  _ dt'  ; 

'  /  t'-t  J  7  f-t 


and  using  nieorem  V  of  Table  1,  we  have 

v2v*(,)  =  1  ^4^  . 


(3.1.4) 


Next  we  multiply  both  sides  of  (3.1.4)  by  1  arid  then  add  (3.1.2)  and  (3.1.4)  ; 


we  then  obtain 


^  ^  v(t)  -  ^ 

c  8  r 


(3.1.5) 


where  V(t)  =  V^(t)  +  1  V^(t)  is  the  anatytic  signal  associated  with  V**!!). 
Thus,  not  only  the  real  disturbance  bit  also  its  associated  analytic  signal 
satisfies  the  wave  equat’on. 

We  recall  that  the  mutual  coherence  function  is  defined  as 


[Toil')  -  <  V,(t  +  T)  V,  *(t)  > 


(3.1.6) 


Differentiating  both  s.i<tes  of  (3. 1.6)  with  respect  to  and  P2  separately, 
formally  interchanging  the  order  of  differentiation  and  integration  on  the  right 
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hand  side,  using  (3. 1.5)  and  tl’.e  stationarlty  of  y(t),  we  obtain  the  two  wave 
equations  (3.1.1). 

3.2  Polychromatic  Fields 

In  this  section  we  shall  study  tlie  limiting  forms  of  the  mutual  coherence 
function  without  making  any  approximations  as  to  the  spectral  width  of  tiie 

illumination.  It  will  be  shown  that  an  optical  field  is  coherent  if  and  only  if  it 
is  monochromatic.  *  Further  It  will  be  shown  that  a  coherent  .field  is  com¬ 
pletely  described,  including  its  coherence  properties,  by  a  simple  wave  function 
(dependlrig  on  one  point  only ). 

It  will  also  be  shown  hi  this  section  that  an  incoherent  field  cannm  exist 
in  free  space  even  the  Illumination  is  polychromatic.  It  is  possible,  however, 
to  define  an  incoherent  aourca  5ih1  the  definition  will  be  seen  to  be  consistent 
with  the  above  consideration  namely,  that  an  incoherent  source,  as  defined  here, 
always  gives  rise  to  a  partially  coherent  field. 

3  2.1  The  Form  of  tor  a  Coherent  Field 

Coherence  la  characterised  by 


yi2(^|  =  1  , 


which  implies  (c.f.  Chapter  1  (1.2)  ) 


Ai  Aj 


e 


(3.2.1) 


(3.2.2) 


In  much  of  the  current  literature  the  term  monochromatic  is 
erroneously  used  to  descrllw  light  with  a  small  but  finite  spectral 
width.  However,  the  term  is  used  here  in  its  strict  sense. 
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Here  we  have  set 

K  -  V  Qs«»  =  1.2)  .  (2-2-3) 


and  (ft  j2(‘^)  is  a  rea^  function  of  tlie  time  delay  t,  and  the  coordinates  of 

_  i  0  i  2  ^  ^ 

Pj  and  P^.  Since  !T2('^)  i®  analytic  signal,  the  function  e 

is  a  unimodular  analytic  signal.  Hence,  according  to  Theorem  VIII  (Chapter 

2)  it  must  be  of  the  form 


r  i'-\ 
12^-^ 


Aj  A2 


i()3  +  2  JT  i’qT) 


'T 

“n’  ■  M 

i 

n=l 

1  j 

“n  1  %  -  M 

.  (3.2.4) 

Here  the  a  are  corr.ylex  constants  witii  complex  conjugates  a  *  ;  and  the 

Imaginary  part  of  each  a^  is  greater  than  zero.  In  (3.2.4)  e 

is,  when  considered  as  a  function  of  a  complex  variable,  z,  a  raeromorphic 

function;  and  the  product  is  taken  c  -er  ali  the  poles  z  =  a  *  .  The  a 

n  n 

are  all  fln^‘e  and  non-zero.  The  constants  0  and  i'q  are  real.  While 
(3.2.4)  represents  the  most  general  unimodular  analytic  slgiial,  it  can,  as 
will  now  be  shown,  ije  interpreted  as  a  niutvai  coherence  function  only  in  the 
degenerate  case  where 


=  Aj  A2  e  ^  .  (3.2.5) 

By  definition  represents  the  complex  cross  correlation  of  the 

disturbances  at  two  points,  and  P2 ,  in  the  field.  Kence,  when  Pj  coin¬ 
cides  with  P2,  the  corresponding  mutual  coherence  function, 
complex  auto-correlation  function.  The  real  part  of 
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real  auto-correlation  function 


[^1(7^)  =  2  =  2  <v/(t+T)  v/(t)>,  (3.2.6) 

(c.f.  Chapter  2,  section  2).  Therefore,  using  also  the  fact  that  the  field  is 
stationary, 

fTl'’'‘(^)  =  »  (3-2-7) 


(c.f.  Chapter  2,  section  2).  This  condition,  (3.2.7),  may  be  expressed  by  the 
stateipent  that 


(t)  sin  2  jTi/T  dr  a  0 


(3.2.8) 


for  all  v.  In  .Appendix  2  of  this  thesis  it  is  shown  that  (3.2.8)  can  be  satis¬ 
fied  if  and  only  if  (3.2.4)  assumes  the  degenerate  form  given  in  (3.2.5),  i.e.  if 


2  -2jri  (^0^  +  ^  ) 
Aj  e 


(3.2.9) 


The  physical  significance  of  (3.2.9)  is  that  all  the  energy  in  the  field 
at  Pj  is  contained  in  the  single  spectral  component  Vq  .  Since  (3.2.9)  is 
valid  for  all  Pj,  it  is  clear  that  the  field  is  everywhere  monochromatic. 
Further,  since  (3.2.9)  is  deduced  solely  from  the  time  dependence,  the 
"constant"  Vq  could  in  principle  depend  on  the  coordinates  of  Pj)  but  ihe 
cross-correlation  of  two  monochromatic  disturbances  of  different  frequencies 
is  ldent..cally  zero,  and  therefore  s/q  must,  in  fact,  be  a  true  constant  charac¬ 
teristic  of  the  entire  field.  Since  the  field  Is  everywhere  monochromatic  of 
frequency  Vq,  the  t  dependence  of  172^^)  must  be  of  the  form 

2iriVnT 


^22  ® 


(3.2,10) 


t 


I 
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where  Fj2  is  a  function  of  the  coordinates  of  Pj,  and  P2  only.  From  the  re¬ 
quirement  that  the  field  be  coherent,  (3.2. 1),  and  the  fact  that  the  field  Ls 
monochromatic,  (3.2.10),  it  follows  that 

I0j2  +  2irt/QT) 


(^2(t)  =  -^1  ^2 


(3.2.11) 


where  i®  ®  function  of  the  coordinates  of  Pj  and  P2  • 

The  form  of  the  function  j3j2  ^  deduced  from  the  fact  that 
satisfies  the  two  wave  equations,  (3.1).  Substituting  (3.2. 11)  into  (3.4)  we  obtain 

[Vg*  +  ]  Aj  A2  =  0  (8=1,2),  (3.2.12) 

where  k  =  2v  t^c.  Performing  the  operations  Indicated  in  (3.2.12),  we  find  that 


w* 

+  k 


+ 


‘ ^12 


V,A 


^1' 


0; 


. (3.2.13) 


and  a  similar  equation  involving  V2  and  A2  holds.  Equating  the  real  and 
Imagink.y  part^  of  (3.2.13)  separately  to  aero,  we  obtain  the  two  equations 


2(V,^J,).  (V,A;) 


(3.2.14) 


and 


(3.2.1:) 


The  right -h&:Ki  side  of  (3.2.15)  is  a  function  of  the  coordinates  of  Pj  only 
and  we  can  therefore  write  (3.2.15)  in  the  form 


f  (Xj,  >J,  Zj)  ■ 


3= 


(3.2.18) 
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On  the  left  side  of  (3.2.16)  the  coordinates  ut  F2  (contained  in  ^^2)  be 
regarded  as  parameters.  Equation  (3.2.16)  it  therefore,  of  the  form  of  the 
eikor.al  equation  of  geometrical  optics  or  the  hamiltcr- Jacobi  equation  of  dy¬ 
namics.  The  general  solution  of  this  equati  n  is  well  known  (c.f.  Born  and 
Wolf  (1959),  page  721),  Let  6®  be  the  value  of  ®  surface  over  which 

U  is  constant  (independent  of  Xj,  y^,  Then  the  solution  is 


(3.2.17) 


^12  =  (p,q,  X2,y2,Z2‘^+  /  ^(x.y.z)ds; 


where  the  integral  is  taken  along  the  extremal  of  the  variational  problem 


6  /  Js  =  0  ; 


(3.2.18) 


0  0 
and  Pj  is  a  typical  point  on  the  surface  '  f'^^ther  p  and  q  are 

two  free  parameters  wUicI  may  be  regarded  as  characterizing  the  orientation 

of  the  surface 12  - 

It  IS  clear  ircni  ’nation  (3.2,17)  that  0^2  b*’  expressed  as 


(3.2.19) 


where  and  ^2  depend  only  (  n  the  coordinates  of  and  P2 ,  respectively. 


Hen.  e  (3.2,11)  becomes 


Aj  e 


^^2  1 


(3.2.20) 


Interchanging  the  roles  of  Pj  and  P2  in  (3.2,20),  and  using  Lemma  I 
(Chapter  2),  we  find 

Aj(Pj)  =  A2(P2)  ,  (3-2.21) 

and 

/3j(Pi)  =  -^(Pj)  +  2n  57  ,  (3.2.22) 

where  n  is  any  integer.  If  we  now  introduce  the  function 

U(P)  =  A(P)e^'^^^^  (3.2.23) 

we  may  rewrite  (3.2.20)  in  the  form 

f;'2(T)  =  U(Pj)  U  (P2)e  ^  .  (3.2.24) 

Thus  in  a  coherent  field  [^2(71  Is  of  the  form  given  by  (3.2.24).  That 
tlj3  converse  is  true,  namely  that  a  mutual  coherence  function  of  the  form  given 
by  (3.2.a4)  always  characterizes  a  coherent  field,  may  be  seen  by  substituting 
(3.2.24)  into  the  definition  of  f72(^'  (1.2.10). 

Thus  with  no  approximation  on  the  spectral  width  of  the  illumination  the 
following  theorems  have  been  established  : 

Theorem  K  :  An  optical  fit  Id  is  coherent  if  and  only  if  it  is  mcnc- 

chromatic. 

Theorem  X:  The  mutiai  coherence  function  for  a  coherent  optical 

field  can  be  expressed  in  the  form  given  by  (3.2.24),  i.e.,  as  a  simple  periodic 

-2ttVQT 

factor,  e  ,  multiplying  the  product  of  a  wave  function,  U,  evaluated 

* 

at  Pj  with  its  complex  conjugate,  U  ,  evaluated  at  Pj  • 


3.2.2  Form  of  (12^'’')  Incoherent  Field 
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Incoherence  is  characterized  by 

Yl2^r)  1^0,  (3.2.25) 

which  implies  (since  the  intensity,  d(®)>  assumed  to  be  finite) 

r^2(T)  3  0  .  (3.2.26) 

However,  by  definition 

f7j(T)  =  I(Pj)  -  <  Vji>-T)  v/(t)>.  (3.2.27) 

TTie  auto-correlation  function,  called  the  self  intensity;  and  it 

is  clear  from  (3.2.27)  that,  apart  from  the  trlvia.l  case  of  an  identically  zero 
field,  the  self  intensity  is  not  identically  zero.  From  (3.2.26)  and  (3.2.27) 
it  follows  that  f  e  mutual  coherence  function  in  an  incoherent  iieid  should  be 
of  the  form 

[72(t)  = 

By  a  simple  generalization  of  the  argument  used  by  the  present  author  for 
quasl-mo.TOChromatlc  illumination*  (see  Parrent  (1959b))  It  will  be  shown  that 

f^2('^)  cannot  be  of  the  form  given  by  (3.2.28)  and  at  the  same  time  satisfy 
the  wave  equation  (3.1.1).  This  theorem  is  most  easily  demonstrated  In  thj 
friKiuency  ctomain.  Accordingly,  we  Introduce  the  mutual  spectral  density, 
f^2(‘')  .  which  may  be  defined  as  the  Fourier  transform  of  ®  ' 

*  Actually  the  argument  used  In  that  paper  Is  slightly  different  from  the  argu 
ment  given  here  but  the  general  structure  of  the  development  is  the  same. 
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(Taf*') 


2ff  tVT 

(72(^6  dT,. 


(3.2.29) 


It  has  beon  shown  (Chapter  2)  tliat,  since  analytic  signal, 

[”*^2(1*)  ■  0  for  negative  v.  By  the  Fourier  inversion  theorem  it  then 


follows  that 


00 

/fe 


-2v  ivT 

(v)e  dv. 


(3.2.30) 


Substituting  from  (3.2.30)  into  (3.1)  we  obtain 


[  *  {^n  iTa*^) 


0  (8  =  1,2). 


In  terms  of  (3.2.28)  can  be  expressed  rs 


f7a(«') 


Pi  >'  Pa 


1  (?,»)  p 


/ 

1  =  P2 


(3.2.31) 


(3.2.32) 


where  .(P,v)  is  the  Fourier  transform  of  I(P,  t). 

Let  V  be  a  finite  volume  of  apace  throughout  which  the  field  ic 
assumed  to  be  incoherent  (i.e.  (3.2.32)  satisfied  throughout  V  )•  Let  .S 
be  any  closed  surface  contained  in  V  .  By  repeated  aj^llcction  of  Green's 
theorem  (see  Parrent  (1959b),  or  Chapter  4)  we  obtain  the  formal  solution  to 
the  Helmholtz  equations,  (3.2.31),  as 


.  (3.1.33) 


where  Gj  and  G2  are  Green's  functions 
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satisfying  (3.2.31)  and  vanishing  over  S;  Sj  and  S2  are  points  on 
Z>  ;  and  and  P2  are  points  within  £  . 

Since  in  (3.2,33)  we  allow  S,  and  S,.  to  explore  the  surface  £  inde- 
pendently,  tiie  integral  is  four-dimensional.  Tlie  integrand  is,  how'ever,  only 

two-dimensiona'  n',^nce  it  is  non-zero  cnly  when  ti'*e  two  points  Sj^  and  S2 

fA 

coincide.  Iheiefore,  if  i  (Sj,S2,  is  assumed  to  he  everywhere  finite  the 
integral  is  ident.  nil,'  ..oto,  i.e  , 


for  ail  Pj  and  P2  (including;  =  p2)-  conclusion,  however,  con¬ 
tradicts  the  assumption  that  (3.2.32)  m  satisiied  throughout  V  We  thus 

find, 

Tlieorem  XI  ;  An  incoherent  field  cannot  exist  in  free  space. 

However,  following  Blanc -Laplerre  and  Dumontet  (1954),  we  may  define 
an  incoherent  source  as  one  for  which  Uie  luut'ual  coherence  functiu.n  is  of  the 
form 

f^2(r)  -  1(P2,t)6{P2  -  Pj)  (3.2.35) 

for  all  pairs  of  points  on  the  source,  where  0  is  the  Dirac  delta  function. 

3. 3  Quasi- Monchromatic  Fields 

The  theorems  derived  in  the  previous  section  are  valid  for  illumination 
of  arbitrary  spectral  width.  However,  in  most  of  Uie  current  applicalions  of 
coherence  theory  one  deals  with  quasi-mcnochrornatic  illumination.  Further, 
as  pointed  out  in  Chapter  1,  most  of  the  results  of  coherence  theory  established 
in  the  literature  are  applicable  only  to  this  special  case.  Therefore,  to  relate 
the  results  obtained  in  this  thesis  to  application  and  to  investigations  reported 
in  the  literature,  it  will  prove  useful  to  examine  in  some  detail  the  limiting  forms 

of  some  of  our  formulae  for  quasl-monochromatic  light. 


V 
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An  optical  field  is  said  to  be  quasi-monocluomatic  if  the  spectral  width, 
A  1/  ,  of  the  light  is  small  compared  to  the  mean  frequency,  F  i.e if 

4^-  «1.  t3.3.1) 

_  r 

More  precisely,  u  and  A ;  ,  may  be  defined  In  terms  of  tlm  mutual  spec¬ 
tral  density  as* 


^  I  A  i  ^ 

J/  I  I 


(3.3.2) 


and 


cc 


(3  .3.3) 


It  can  be  shown  that  when  (3.3.1)  Is  s.ati><fied  the  mutual  co.hererce  func 
tion  may  be  expressed  in  the  form  (o.f.  Appendix  2) 


-2  .flrr 


(1  tI< 


A  ^ 


(3  3,4) 


provided  we  restrict  our  attention  to  sufficiently  email  tirro  dlfierence  j,  s  in 
dicated  on  the  right  in  (3.3.4).  Quasi-nionochi\)matlc  fields  may  he 
*0  be  coherent  if  the  condition 


yi2(^ 


(1 


T  <.v: 


-v) 


(3.3.5) 


is  satisfied.  Adoption  of  this  considerably  weaker  condition  is  cquicaient  to 
♦  Following  Wolf  fl958b). 


<4 
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l-r- 

■  s'*  • 


examining  on  v  th>?  central  iringes  in  a  Young'.?  interference  experiment  used 
to  measure  the  degree  of  coherence,  ‘^ich  a  "modified  dafinitioa"  is  reason¬ 
able  for  quasi  monochromatic  fieids  since  monochromatic  ilburdnation  is  an 
unrealizable  idealizauon.  and  light  of  sufflctently  narrow  spectra!  width  behaves 
ab  coherent  light  under  suitaMe  conditions.  ♦ 

Substituting  from  (3.3. -i,  into  the  wave  eq'iation''  (3.1.1),  we  obtain 


Vg“  f  gH  --  0.  (s  =  l,2) 


(3.3.6) 


where  k  i.s  the  mean  wave  number.  Combining  (3.3.4j  and  (3.3.5)  we  find 
the  mui..:ai  inienslty  for  coherent  uuasi-monochronatic  fields  to  be  of  the  form 


»*!!►-  ! 


C." 


Toio) 


Aj 


(3.3.7) 


w.h-.: re  s  before  -  A„.  f  rom  (3.3.6)  and  (3.3.7)  and  the 

f 

argumentb  of  Section  ?.  of  thrs  chapter  we  obtain  the  form  of  Theorems  DC  and 
X  for  the  quasi-ujonochromatic  approximation 

Theorem  XII  .  A  quasl-monochronnuic  field  is  coherent  if  and  only  If 

the  mutual  lnter»slt>'  can  be  expressed  as  ihe  product  of  a  wave  function,  U, 

* 

evaluated  at  with  us  coiiiplex  conjugate,  U  ,  evaluated  at  P2,  i.e.,  if 


M9^) 


U(Pj)U  (P2) 


^3.3.8) 


I'i' 


U(P)  is  the  function  defined  In  Section  2. 

This  result  is  in  igreenreiit  with  the  observation  that,  (or  the  purpose  of 
cah  u’ntlng  the  intensity  In  diffraction  phenomena  involyi;tg  short  path  differences 


We  ix)bit  out  that  quasi- n'iOi>ocnroinaiic  light  may  also  be  incoherent  or 
partially  coherent,  c.f.  Chapter  4  of  this  thesis;  on  the  other  hand, 
monochromatic  light  is  always  coherent. 


r 
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one  may  tr  eas  coherent  quaai-tnonochromatic  fields  as  monochromatic  That  la 

a  coherent  quaei- monochromatic  field  completely  specified  ly  a  single  mono-  ^  j 

chri)matic  wave  function  . 

The  tneorem  of  the  preceding  section,  which  states  that  an  incoherent 
optical  field  cannot  exist  in  free  space,  is  valid  in  general  and  hence  is  obviously  | 

true  for  quasi-monochromatic  fields,  'ihe  theorem  may,  however,  be  demonstrated  I 

I 

directly  for  the  case  considered  here  (this  reoult  was  originally  obtained  for 
quasi- monochromatic  fields  in  this  manner  (see  Parrent  (1959b)).  An  incoherent 
quasi' monochromatic  source  may  be  defined,  by  analogy  with  the  general  case, 
as  a  source  the  mutual  intensity  of  which  is  of  the  form 

[^2(0)  =  1(82)6(82  '  Sj)  .  (3.3.9) 


I 

I 
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CHAPTER  4 

THE  PROPAGATION  OF  PARTIALLY  COHERENT  LIGHT 

The  laws  governing  the  propagation  of  partially  coherent  opti'^al  fields 
have  been  of  central  importance  in  the  development  of  coherence  tlieory.  The 
most  important  of  the  earlier  contributions  to  this  aspect  of  the  theory  is  a 
theorem  discovered  separately  by  van  Cittert  (1934,  1939)  and  Zernike  (1938). 
Under  suitable  conditions  on  the  spe'trai  width  of  the  illumination  (quasi- 
monochromatic  light),  this  theorem  expresses  the  mutual  intensity  on  a  plane 
illuminated  by  an  incoherent  plane  source  in  terms  of  the  intensity  distribution 
across  th  source.  In  the  same  paper  Zernilco  also  derived  an  approximate 
law  for  the  propagation  of  the  mutual  intensity.  Later  Hopkins  (1951)  derived 
thes^  theorems  in  a  different  way. 

Li  the  formulation  of  coherence  theory  used  here  these  tlieorems  will  b*' 
seen  to  be  limiting  or  aporoxlmatc  forms  of  a  Green's  function  solution  to  the 
wave  equations  which  were  derived  in  Chapter  3.  hi  the  present  ctiapter  we 
will  OihJ  .n  Ihe  general  solution  for  the  propagation  of  mutual  cohcroncc  from 
a  plane  polychromatic  partially  coherent  source.  The  limiting  for.ms  of  this 
general  solution  will  be  examined  in  some  detail  both  for  sources  of  wide  and 
small  spectral  ra:igcs,  and  tlie  van  Cittert -Zernike  theorem  will  be  shown  to 
represent  an  approximate  form  of  the  incohc’-ont  limit  of  the  quasi-monf^chromatic 
solution. 

4.1  General  Solution  for  thejpistribution  of  MvHoal  QQherpr,ce  froril  a  Pl^l? 

Source 


» 


In  tills  .seetJer  we  shall  detcrmiiu'  tiie  mutual  coheionce  function  for  a 
field  created  by  a  plane  polychroniatJ'.'  source,  hi  Figure  2,  which  serves  to 
define  the  cov>rd'nates,  a  is  the  plane  containing  ihe  extended  polychromatic 


source  with  a  known  distribution  of  mutual  coherence.  P,  and  P2  are  points 
In  the  illuminated  field  and  S.  and  are  points  in  the  plane  of  the  Kjurce. 

I  t 

Pj  and  Pg  are  the  mirror  images  of  the  points  P^  and  Pg  respectively  in 
the  plane  7  . 

As  shown  in  Chapter  3  the  propagation  of  the  mutual  coherence  function 
In  vacuum  is  governed  by  the  two  wave  equations 


Figure  2 

We  assume  that  rj2(^)  is  known  for  ail  pairs  of  points  Gj  and  S2  in  the 
plane  a  Let  ^  Fourier  time  transform  of  !72(  '')  •  Since 

is  an  analytic  signal  (see  Chapter  2),  li  contains  only  positive  fre- 
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I 


-2ir  ivT 


di^. 


(4.1.2) 


anc:  by  the  inversion  theorem 


(lafr) 


«0 

J"f72(T)e 


2  ttIvt 


0 


dr  V  >  0, 


V  <  0. 


(4.1,3) 


Substituting  from  (4.1.2)  into  (4.1.1)  and  interchanging  the  order  of  lists- 
gration  and  differentiation  we  obtain 


«0 

J  +  k^(  u)  ]  |52(»')e 


-2)r  ivT 


dv  =  0  (s  =  1,2) 


(4.1.4) 


since  (4.1.4)  must  hold  for  all  r,  we  have 


1  f72(‘^)  =  (s  =  1,2) 


(4.1.5) 


where  k(v)  =  2r  v/c.  Thus,  ea  h  spectral  component  of  satlsflss 

the  two  scalar  Helmholtx  equations,  (4.1.5). 

Equation  (4.1.6)  can  be  formally  solved  by  employing  Green's  fusctlons. 


To  this  end  we  integrate  first  over  the  coordinates  of  and  obtain* 


p  r  ^ 

np,,a,,  >•)  =  -  S  /  ^  (\h’ 


»G, 


V)  dSj  . 


(4.1.6) 


•  (^2(7)  aiid  f5|'2(v)  will  be  written  as  nPi>P2»^^  ^ 

when  necessary  to  stress  the  space  dependence. 


»') 
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Here  Gj  is  a  Green's  function  satisfying  the  equation 

[Vg2  +  k2(t-)]  Gj  0  (4.1.7) 

with  the  boundary  condition 


G(Sj)  =  0. 


(4.1.8) 


Equation  (4.  A. 6)  provides  the  boundary  condition  for  the  solution  of  the  second 
Helmholtz  eqxiation;  and  employing  the  same  theorem  again  we  obtain 

l^{Pj,Pj,^)=-  r(P,,S2,^)  dSj,  (4.1.9) 

a 

by  integrating  over  the  coordinates  of  S2,  where  G2  is  a  second  Green's  func¬ 
tion  satisfying  the  same  conditions  .as  G^  .  Substituting  from  (4.1.6)  into 
(4.1.9)  we  obtain 


A 


(P^Pj,!-) 


cr  a 


dS^  dS2- 
(4.1.10) 


In  order  to  determine  the  exact  form  of  Gj  and  Gg  for  cases  of  physi- 

A 

cal  interest,  we  impose  on  radiation  condition  of  Sommerfeld. 

In  essence  this  condition  implies  that  the  distant  field  is  essentially  that  of  a 

diverging  spherical  wave.  By  applying  the  radiation  condition  to  Vj(j^)  and 
♦  A 

Vn  (i^)  and  appealing  to  the  definition  of  [72^^)  terms  of  these  functions,  it 

A 

is  readily  seen  that  must  behave  asymptotically,  lorki-j  >>  1  and 

kr.j  >  >  1  as 

£1 


172  (i'')  l'^2^ - 

*’l‘’2 


(4.1.11) 
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where  9^  and  are  the  spherical  polar  coordinates  of  Pg.  To  satisfy 
this  condition  and  those  defining  Gj  and  G2  the  required  Green's  func¬ 


tions  are 


(4.1.12) 


G„  =  5 - ^  (4  113) 

*  r  r 

f  t 

where  ri,  r2,  ri  and  rg  are  defined  in  Figure  2.  That  these  are 
the  required  Green's  functions  can  be  seen  by  direct  substitution.  Before 
substituting  Into  (4.1.10)  we  obtain  the  normal  derivatives;  we  have 


(ikrj-l)  f 


^^1  8r, 
e  1 


—  +  (I'ikr, ) 

n  '  1 ' 


'  s  e. 


r,'  8n 


(4.1.14) 


Noting  that  r .  j 
^  icr 


r.  and  that 
^  o 


0n  a 


and  setting  cos  6  -  z  /r  we  may  rewrite  (4.1.14)  as 

S  o  S 


0  G, 

^  =  2(lkr,-l) 


COS0,  — 


and  similarly 


-  -  2(lkr  +  1)  008^2’ 


(S. i  io) 


(4.1.16) 
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On  substituting  from  (4.1.15)  and  (4.1.16)  into  (4.1.10),  we  obtain 


/hsj.Sj,!/)  (1-Uirj)a+ 

P')  a 


lk(rj-r2) 

ikr^).  cos  £L  cos  $2  - -  dS2 


.  (4.1.17) 

Equation  (4.1.17)  is  the  contribution  from  a  single  spectral  component, 

V.  The  complete  eolutlon  is  obtained  from  (4.1.17)  by  taking  the  Fourier 
transform  of  both  sides  of  (4.1.17),  i.e. 

«o  llc(r  ~T  ) 

^2^^  yyji^8i.S2,»')  (l-lkrj)(l+ikr2)post?jCos«2^ -  ®  dSjdSgdi' 


0  o  a 


>•1^2 


. (4.1.18) 

Since  dSj,  dSj,  and  d^  are  independent,  we  may  invert  the  order  of  integra¬ 
tion  and  obtain 


GjW  .  -Lj.  ff  cos  0 


where 


cos®  2  0  j2(^  dSj  dflj  , 

.  (4.1.19) 


,  -2liu(  T  -  — ^ 

O  lj(^  "  Jr(Si,S2,y)e 


r.  -  r. 


-) 


dp 


0 

so 


r,  -r. 


.‘(Xli)  f  2,c  T/S  o 


9 

«0 


(Si,S2,p)e 


dp 


2  2  A  .  ,  -2»ip(T  -  -1- 
^  dp 
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We  recall  the  following  well-known  theorem  from  Fourier  analysis.  If 


fit)  = 


r 


-27rit't 

g(i/)e  di^ 


(4.1.20) 


then 


=  f  (-2irii/)"  g(;^)e 

3t  J 


-2jrii't 


d^■ 


(4.1.21) 


Using  (4.1.21)  may  be  evaluated  giving 


p  /  ^  ^  ^ _ —  ^1^2  A2^/SitSniT_  1  2 

ni2(r)  =r  [Si,S2.r  -  P  ^  c  ^  P  2  _ 


.  .  .  .  (4.1.22) 


and  the  final  solution  la 


^12^’''  -  jj  cos  cos  ^2  ^^r(Sj,S2.  T  - 


cr  cr 


dSj, dS2 , 

.  .  (4.1.23) 


where  ^  Is  the  differential  operator 


1  + 


^1  ■  *'2  3 

5  Tt 


■■1^2  8^ 


(4.1.24) 


Equation  (4.1.23)  is  the  general  solution  for  the  mutual  coherence  in  a 
field  produced  by  a  plane  polychromatic  source. 

4.2  Coherent  and  Incoherent  Polychromatic  Sources 

In  this  section  we  examine  tlie  limiting  forms  of  (4.1.23)  for  colyerent 
and  incoherent  sources.  No  approximation  on  the  spectral  width  of  the  illumina¬ 
tion  will  be  made  here.  We  will  show  that  1.)  a  coherent  source  always  gives 


m 


\ 


rise  to  a  coherent  field  and  2.)  an  incoherent  source  always  creates  a  partially 
cohere.it  field.  By  examining  the  incoherent  limit  we  will  obtain  a  rigorous 
generalization  of  the  van  Cittert-Zernike  theorem  to  polychromatic  sources. 
4.2.1  Coherent  Source 

It  was  shown  in  Chapter  3  that  coherence  is  characterized  by  a  mutual  co¬ 
herence  function  of  the  form 

p  *  -2niVf.T 

I  (SpSg,  T)  ^  U^(Sj)  (82)6  ^  .  (4.2.1) 

Taking  tlie  Fourier  transform  of  both  sides  of  (4.2.1)  we  obtain 

[^(8^,82,!^)  =  U^(Sj)  U*(82) 

Substituting  from  (4.2.2)  into  (4.1.17)  we  obtain 

A 

r(Pj,P2,i^)  =  U(Pj)  U*(P2)5(i^-V  ’  (^•2-3) 

where 

U(Pj)  = 

a 

Taking  the  inverse  transform  of  both  sides  of  (4.2.3),  we  obtain 

n  ♦  -2jrii/„T 

I  (Pj,P2,t)  =  u(Pj)U(P2)e  (4.2.5) 

From  (4.2.5)  and  Theorem  IX,  p.  57,  it  is  clear  that  the  following  theorem 
holds  : 

"nieorem  XIII  ;  In  vacuum  a  coherent  source  will  always  give  rise  to 


(Sj)  cos  1?  (1  -  ikrj)  £_ 

r  i 


ikr, 


dS, 


(4.2.4) 


a  coherent  field. 
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4.2,2  Incoherent  Source 

In  Chapter  3  it  was  shown  that  an  incoherent  field  cannot  exist  in  free 
space  although  an  incoherent  soi:rce  may  be  defined.  In  the  present  section 
we  will  show  that  the  definition  of  such  a  source  is  consistent  with  the  result 
that  an  incoherent  source  will  give  rise  to  a  partially  coherent  field. 

By  definition,  an  incoherent  source  is  characterized  by  a  mutual  coherence 
function  of  the  form 


1^.  , 
§■'.  ’  '  * 


r(Sj,S2,T)  =  I(S2,t)5  (Sg  -  Sj) 

A 

Hence  f*  (8^,82,  v)  is  given  bj' 


A  ^ 

r(Sj,s  v)  =  i(Sj,v)a  (Sj  -  Sj). 


(4.2.6) 


(4.2.7) 


I. 


We  substitute  (4.2.7)  into  (4.1,17);  and,  after  integrating  over  82,  we 


obtain 


’pPa,*-)  4i(s,, 


v)  (1-ikrj)  (l+ikrg)  co8  0j  cos^g  s. _ 


lk(rj-r2) 


. (4.2.8) 


where  r^  and  r2  are  now  interpreted  as  the  distance  from  S  to  and  P2 
respnctivcly. 

By  substituting  from  (4.2.8;  into  (4. ].x8),  and  Integrating  over  Sj,  we 
obtain  the  mutual  coherence  function  for  the  field  of  an  extended.  Incoherent, 


polychromatic,  source  , 


^12^^  =  J'coae^  cos  0  2 


"l  -  *-2 


)]  dS  . 


. (4.2.9) 


I 
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Equation  (4  2.9)  is  the  generalization  of  the  van  Cittert-Zernike  theorem 
for  polychromatic  sources.  It  expresses  the  mutual  coherence  function  under 
the  coiiditions  stated,  in  terms  of  the  self  coherence  across  the  source. 

Since  an  optical  field  is  coherent  if  and  only  if  It  is  monochromatic,  it  is 
cleaii  that,  while  an  incoherent  source  gives  rise,  by  (4.2.9),  to  a  partially  co¬ 
herent  field,  an  incoherent  source  cannot  create  a  coherent  field. 

4.3  Quasi- Ktonochromatic  Sources 

As  mentioned  earlier,  most  of  the  current  applications  of  coherence  ttu^ory 
involve  quasl-monochromatic  light.  For  this  reason  we  shall  treat  the  quasi- 
monochromatic  limits  as  separate  problems  deducing  them  directly  from  the 
general  solution,  (4.1.23),  rather  than  obtaining  them  as  special  catses  of  the 
results  of  the  previous  section. 

4.3.1  Coherent  Source 

Since  in  the  quasi-monociiromatic  approximation  a  field  is  described  as 
coherent  if  the  condition  |  i  ^  satisfied  only  for  sufficiently  small 

T,  coherence  Is  a  considetably  weaker  condition  for  the  class  of  problems  ('on- 
sidered  tha.:  in  the  general  case. 

A  coherent  quasi-monochromatic  source  has  the  property  that  for 
[  T  I  "^<2^  ifs  mutual  coherence  function  is  of  the  form 


r(Sj,S2,i  =  U^(Sj)  .  (4.3.1) 

3u''8tituting  from  (4.3.1)  Into  (4.1.23)  yields 

r(Pi,P2,t)  =  U(P,)  U*(P3)e'^'^*'’‘  (|rl«  £)  ,  (4.3.2) 


where 


U<Pi) 


J 


u 


(Sl)(l  - 


ikr. 


ikrj)cosflj  ^ 
r 


dS 


(4.3.3) 


1 


57 


Equations  (4.3.2)  and  (4.3.3)  constitute  the  mathematical  restatement  of 
Theorem  XIII  for  quaai-monochromatic  fields. 

Tliis  result  is  not  simply  a  s^iecial  case  of  the  previous  result  even  though 
in  Section  4.2  we  were  able  to  establloh  the  tlieorem  with  no  apprcriniations.  rJ'e 
result  is  perhaps  more  interesting  and  certainly  of  more  immediate  practical 
importance  in  the  present  case.  In  the  general  treatment  where  no  approximations 
on  the  spectral  width  of  the  light  were  made  we  concluded  that  an  optical  field 
is  coherent  if  and  only  if  it  is  monochromatic.  It  is  certainly  to  be  expected 
that  such  a  field  will  remain  coherent  as  it  p^  cnagates  in  free  space.  It  is  also 
to  be  expected  that  such  a  field  can  be  completely  specified  (apart  from  polariza- 
iion  effects)  by  a  simple  wave  function,  depending  on  the  coordinates  of  one 
point  only. 

The  fact  that  these  results  hold  under  the  narrow  spectral  width  approxi¬ 
mation,  however,  indicates  something  more.  It  is  evident  from  (4.3.2)  and 
(4.3.3)  (and  in  fact  from  the  results  of  Chapter  3),  that  if  we  restrict  ourselves 
to  phenomena  involving  sufficiently  smal?  patli  differences,  (It I  ® 

quasi-  tonochromatic  field  may  behave  in  some  respects  like  a  monochromatic  field. 
But  there  is  an  essential  difference  between  tnese  two  types  of  illumination  : 

A  monochromatic  field  is  evorywhure  coherent  for  all  t,  while  a  quasi- 
monochromatic  field  cannot  be  coherent  In  the  strict  sense,  j  j  "I* 

but  only  in  terms  of  the  weaker  condition,  |  I  *  ^  • 

Further  a  quasi-monochromatic  field  may,  as  will  be  discussed  in  Section  4.3.2, 
be  incoherent. 

In  spite  of  these  considerations  the  similarity  between  these  two  types  of 
fields  for  a  large  class  of  diffraction  phenomena  has  led  to  a  loose  usage  of  the 
term  monochromatic  and  hence  to  such  meaningless  statements  frequently  found 


r 
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in  the  literature  "partially  coherent  {or  incoherent)  monochromatic  light"  and 
the  phrase  "to  extend  the  concepts  of  partial  coherence  to  light  with  finite 
spectral  width". 

4.3.2  tocQJifirftai  Sonri^ 

Under  the  conditions  justifying  the  quasi-monochromatic  approximation 
we  may  write  (see  Appendix  2) 


(|H  «:^)  •  (4.3.4) 


On  substituting  from  (4.3.4)  into  (4.1.8)  and  integrating  over  u,  we  obtain 


-2x11^7  ff 

- 2 —  I  I  (^"^^1^  +  cos  6,  cos  C.^\  (SpS2,0) 

(2tt)  JJ 


n  <y 


Ikfrj-rj) 


68^682- 


. (4.3.5) 

Equation  (4.3.5)  expresses  the  mutual  coherence  function  for  small  j  t  j  in 
the  field  created  by  a  plane  partially  coherent  qimsi-monochromatlc  source. 

We  now  consider  the  limiting  form  of  (4.3.5)  when  the  source  is  Incoherent. 
By  definition  in  Incoherent  quasi -monochromatic  source  may  be  specified  by  a 
mutual  coherence  function  of  the  form 

rj2(T)  =:  I(S2)5(S2'Sj)e'2^*^^  (h!«^;)- 


- (4.3.6) 


Subrtltuting  (4.3.6)  into  (4.3.5)  and  integrating  over  S2  yields 
■2irl.“r 

- 1  1(8)  (1-lkr,)  (l+ikr,)  cos  0,  cos^o 

2 


lh(ri-r2) 

n9(T)  =  ■® — Tj — /  1(8)  (l-lkr.)  (l+ikr,)  cos  0,  CO8  09  dfl 

.A  j2ir)^  J  ^  1  ^  ^ 


ij.' 


(1.3.7) 


Here  Fj  and  r2  are  now  interpreted  as  the  distance  from  a  typical  point  S  on 
the  source  to  the  field  points  Pj  and  F  respectively,  and  cos  « 

(s  *  1,  2).  Equation  (4.3.7)  expresses  the  mutual  coherence  function,  Ii2(^^  ’ 
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for  sufficiently  small  |  t  |  in  terms  of  tlie  intensity  distribution  across  the 
source. 

If  (4.3.7)  is  evaluated  at  r  =  0,  the  obliquity  factors  are  ignoredand 
attention  is  limited  to  field  points  on.  a  plane  parallel  to  a  we  obt'^in  the  theorem 
due  to  van  Cittert  (1934)  and  Zernike  (1.933)  already  referenced,  viz., 


I 

i 


I(S)- 


ik(rj-r2) 


’'1^2 


dS. 


(4.3.8) 


This  theorem  expresses  the  mutual  intensity  on  the  Illuminated  plane  in  terms 
of  the  intensity  distribution  of  the  incoherert  source. 

in  most  applications  one  is  interested  in  the  form  of  (4.3.8)  in  the  Fraunhofer 
approximation.  The  right-hand  side  of  (4.3.8)  then  reduces  to  a  Fourier  trans- 
fo.^m  of  the  Intensity  distribution.*  If  the  intensity  distribution  is  suitably  norma¬ 
lized  (4.3.8)  becomes 


^12(0)  =  JJu^,  ^  .  (4.3.9) 


where 

^  i  A 

P  =  TT,  ■  ?r-  - 

and 

‘-8  =  V*"'  ' 

,  ^1  ^2 
?  T 


These  coordinates  are 


defined  in  Figure  2,  p.  67.  In  the  region  where  (4.3.9)  Is  valid  R,  and  R,  may 
1)6  taken  as  equal  Equation  (4.3.9)  is  the  mo?)t  commonly  used  form  of  the  van 
Cittert-Zernike  theorem.  We  see  ttet  formula  (4.3.5)  <s  a  generaliza¬ 
tion  of  this  theorem  to  partially  coherent  (but  quaaj-monochromatlc)  sources. 


Actually,  the  validity  of  life  Fourier  traniforr  relation,  (4.3.8),  doss  not 
depend  on  the  Fraunhofer  apDrcxlmatlors  lion  (4.3.9)  ly^ies 

ever  the  obliquity  factors  ana  me  variAticnu.  ^  i  ^  in  (4.3.7)  can  be  ig- 
noted  while  the  Fraunhofer  • '.xln.ar.on  in  ackiitlon  to  tt.ese 

conditions  '.hat  R  >>  --  u  a 


C!IAPTER  5 


IMAGING  OF  EXTENDED  POLYCHROMATIC  SOURCES 
AND  GENERALIZED  TRANSFER  FUNCTIONS 

In  thij»  chapter  we  snail  apply  the  theorems  and  results  of  the  earlier 
chapters  to  the  determination  of  the  relation  between  object  and  image  for 
systems  which  image  extended  polychromatic  objects.  We  shall  treat  the 
problem  prlma^'lly  in  the  spatial  frequency  domain,  an  approach  introduced 
by  Duffieux  (1948)  in  1346.  Since  its  introduction,  the  frequency  domain  analy¬ 
sis  haf?  proved  veiy  poweriul  In  the  t-tudy  of  imaging  systems.  Li  this  r nalyais 

the  imaging  system  is  described  by  a  transfer  function  (also  called  modula¬ 
tion  function,  transmission  factor,  transmission  funciion.  contrast  ren.  tlon 
function,  frequency  response  function).  The  imaging  problem  is  then  solved  as 
follows  :  the  object  and  image  are  described  in  terms  of  the  distribution  of  a 
isuitable  physical  characteristic  of  the  optical  disturbance,  which  characteristic 
is  determined  by  the  degree  of  coherence  of  tlie  object  illumination.  For  exampitj 
an  incoherently  Illuminated  object  is  described  in  terms  of  'he  Intensity  distri¬ 
bution  across  it.  The  spullal  spectrum  of  the  image  is  Uien  obtatne  l  as  the 
product  of  the  transfer  function  with  the  spatial  spectrm.i  ot  uie  object,  ^.e., 
the  spatial  Fourier  transform  of  the  above-mentioned  distribution.  That  is, 
the  optical  system  is  treated  as  a  spatial  frequency  filter.  After  its  introdti:- 
tlon  into  the  study  of  optical  systems,  the  .sfer  function  anclysife  was  a^'plled 
to  the  stix'v  of  the  mapping  problems  oi  radio  astronomy  (see  Bracewe*!  and 
Roberts  (1954)  )  and  radai. 

This  analysis  Ic  particularly  promising  in  the  study  of  cascaded  sysu  s 
as  exemplified  by  Schade's  (1948)  treatment  of  television  systems.  In  cascaded 
systems  the  final  Image  in  the  frequency  domain  is  obtained  by  multiplying  Lii? 
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'ipi?ctrum  of  ^he  object  v/ith  the  product  of  the  transfer  fuucUoiiS  describing 
et'ch  stage  of  the  system. 

The  essential  step  in  this  approach  is  the  lecognition  that  many  c'ptical 
systems  may  to  a  good  approximation  be  treated  as  linear  stationai7  systems 
in  terms  of  their  spatial  as  well  as  theii  temporal  dependence.  All  the 
advantages,  familiar  to  electronics  engineers,  of  performing  linear  system 
analysis  in  the  frequency  domain  may  then  be  realized  in  optical  imaging 
problems.  In  optics  ciiti-e  are,  howevei,  3ome  difficulties  confronting  this 
approach. 

For  example,  functions,  describing  objects  anc’  images,  met  in  the 
analysis  of  optical  systems,  depend  on  spatial  as  well  as  temporal  coordinates; 
some  optical  systems  of  practical  interest  are  not  ".stationary"  in  their  spatial 
variation.  'Ihe  most  significant  difficulty  is  the  fact  t.hat  the  form  of  the  trans¬ 
fer  function  is  determined  by  the  degree  of  ccherer.ce  of  the  object  illumination. 
As  we  shall  .see  i.i  the  following  di /elopment,  some  of  these  difficulties  do  not 
appear  in  analogous  problems  in  radio  astronomy  and  radar,  though  these  iields 
present  other  p.  oblems. 

The  realization  that  imaging  systems  can, under  suitable  conditions,  be 
analyzed  as  linear  stationary  systems  suggest.s  stroiigly  the  application  of  the 
techniques  of  information  and  communication  theory.  Here  again,  ho/'f>ver, 
several  basic  difficulties  arc  encountered. 

Apart  from  the  considerations  already  mentioned,  an  imaging  device  is 
not  in  general  a  communication  sy.nten)  (since  no  opportunity  of  encoolng  the 
input  exists)  but  rather  an  observ.'itlo.i  system.  Tne  difficulties  con  routing 
the  analysis  of  such  a  system  in  terms  cf  information  and  communication  Dieory 
are  discussed  and  illustrated  by  Woodward  (1953)  and  lie  outside  the  domain  ol 
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our  present  discussion. 

In  spite  of  these  considerations  some  results  have  been  obtained  from 
the  application  of  information  theory  to  imaging  systems.  The  chief  contri¬ 
butions  of  this  theory  to  the  study  of  image  formation  are:  1.)  the  demon¬ 
stration  that  an  optical  image  has  a  finite  number  of  degrees  of  freedom  (see 
Fellgett  and  Linfoot  (1955)  and  Gabor  (1956)  )  ;  and  2.)  the  demonstration  that 
the  criteria  for  judging  the  quality  of  an  imaging  device  must  take  account  of 
the  objects  that  the  system  is  to  image  (see  Fellgett  and  Linfoot  (1955)  and 
Schade  (1918)  ).  The  ,  of  these  conclusions  follows  Imrcediately  from 
Shamion's  Sampling  Theorem  (c.f.  V/oodward,  1953)  and  the  fact  that  an 
imaging  system  behaves  as  a  low  pass  filter  with  a  finite  cut  off  frequency.  The 
second  consideration  will  become  evident  from  tlte  subsequent  discussion  of  tills 
chapter. 

Some  of  tlie  vr  rJouc  quality  criteria  for  imaging  devices  which  were  intro¬ 
duced  by  Fellgett  anti  Linfoot  and  Schade  have  beer  evaluated  for  aberrated 
optical  systems  (see  O'Neill  (1956),  Fukul  (1957)  and  Parrent  and  Drane  (1956)) 
and  for  ;  .tenna  tiystems  employing  Doiph-Tchebbycheff  apodlzation*  (see  Drane 
(195?)  ).  Tlie  results  obtained  in  each  case  were  in  good  qualitative  agreement 
with  experience. 


It  can  be  show'n  (see  Oolph  (1946)  )  that  if  the  currents  in  the  elements 
of  a  linear  (antenna)  array  are  prooorLlona!  to  the  coefficients  of  the 
Tcbetejycheff  polynomials  the  resulting  diffraction  pattern  has  Uie  mini¬ 
mum  possible  aide-lobe  level  for  a  given  beam  widUi.  These  polynomials 
provide  a  means  of  varying  the  apodlzatlon  continuously  from  edge  illumi¬ 
nation  (which  gives  the  cosine  squared  diffraction  pattern  of  a  simple 
interferometer  and  hence  the  minimum  beam  width  for  a  given  aperture 
ske)  to  a  binomial  distribution  of  currents  which  gives  a  diffraction 
pattern  consisting  of  a  main  lobe  with  no  side  lobes. 
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In  spite  of  the  difficulties  mentioneo,  the  transfer  function  analysis  has 
contributed  to  the  understanding  of  the  problems  of  image  formation.  Among 
the  interesting  consequences  of  the  application  of  the  transfer  function  analysis 
to  optical  systems  are  the  spatial  filtering  techniques  developed  by  O'Neill 
(1966)  and  Marechal  and  Croce  (1963)  for  sharpening  blurred  edges  and  recover¬ 
ing  wanted  detail  from  an  image  containing  noise  (e.g.  photcgraplilc  grain). 

The  computation  of  the  transfer  function  of  any  given  system  often  leads 
to  numerical  integration;  and  for  this  reason  detailed  analysis  has  been  limited  to 
relatively  simple  systems,  utilizing  strictly  coherent  or  strictly  incoherent 
illumination,  with  small  aberrations  (see  Steel  (1953)  or  witl.  a  single  aberration 
^Parrent  (1956),  Hopkins  (1956),  De  (1955),  O'Neill  (1956)  ). 

The  limitation  to  coherent  or  incoherent  illumination  is,  however,  more 
basic  and  stems  from  the  following  considerations  :  l.)  A  system  imaging  an 
incoherently  illuminated  object  may  be  regarded  as  linear  in  Intensity;  2.)  A 
system  imaging  a  coherently  Illuminated  object  may  be  regarded  as  linear  in 
amplitude;  3.)  Systems  using  parUally  coherent  iHumLiation  are  linenr  in  neiUjer 
of  t.ese  quantities.  H.  H.  Hopkins  (1956)  end  Dumontet  (1954)  extended  the  trans¬ 
fer  function  analysis  to  systems  imaging  partially  coherent  objects  by  showing 
that  such  systems  may  be  regarded  as  linear  ir.  mutual  iiitensity.  The  transfer 
functions  for  systems  with  small  aberrations  and  partially  coherent  objects 
have  been  computed  by  Steele  (1957)  using  the  Hqf^lns  formulation. 

In  tjvery  case,  coherent,  psriiaily  coherent  or  incoherent,  the  analysis  has, 
however,  been  limited  to  quaoi-monochromatlc  light.  It  la  our  aim  in  this  chapter 
to  provide  the  framework  for  the  analysis  in  terms  of  the  transfer  function  of 
syatema  employing  polychromatic  illumination;  and  further  using  the  theorems 
and  results  of  the  preceding  chapters  to  verify  that  the  familiar  solutions  to  the 
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quaui-monochromatic  imaging  problems  follow  from  our  general  solutions  as 
appioxlmate  forms.  It  is  hoped  that,  in  addition  to  providing  the  mathematical 
frameworic  for  Ihe  solutions  to  these  general  problems,  this  approach  will  elim¬ 
inate  the  frequently  encountered  confusion  concerning  the  significance  and  in¬ 
terpretation  of  the  various  transfer  functions.  * 

S.i  General  Formulation  of  the  Imaging  Problem 

hi  this  section  It  will  be  shown  that  the  general  scalar  imaging  problem. 
Involving  partially  coherent  polychromatic  objects,  can  be  completely  solved 
In  terms  of  the  observables,  172^^  >  ^  object  and  image  space  with  no 

recourse  to  the  random  disturbance,  V^(t) ,  itself.  By  dealiiig  solely  with 
the  mutual  coherence  function,  and  functions  simply  derivable  from  it,  our  entire 
analysis,  apart  f^om  the  limiting  forms,  will  involve  only  square-integrable 
functions.  The  advantages  of  such  an  approach  are  immediately  obvious  since 
we  deal  extensively  with  both  Fourier  and  Hilbert  transforms  and  their  respec¬ 
tive  inversion  theorems. 

It  will  be  shown  hero  that  using  this  general  solution  il:  is  tc 

define  generalised  transfer  functions  whose  properties  are  formally  similar 
to  those  of  the  quasl-monochromatic  transfer  functions.  Further  it  will  be 
shown  that  these  new  functions  are  simply  derivable  from  the  aperture  illumina¬ 
tion  function  (pupil  function)  of  the  imaging  system. 


This  confusion  is  discussed  and  illustrated  by  F.  J.  Zucker  in  his  summary 
comments  jHibUshed  in  the  Electronics  Research  Directorate,  Air  Force 
Cambridge  Reeearch  Center,  ARDC,  January  1057. 


5.1.1  K  theinatical  Conventions  and  Notations 

Li  1  :e  devekpment  that  follows  extensive  use  is  made  of  multi-dimensional 
Fourier  transforms.  To  prevent  the  equations  from  becomiiig  too  unwieldy  the 
following  conventions  and  condensed  notations  will  be  used. 

Cartesian  coordim  tes  will  be  denoted  by  ( I ,  rj )  in  object  space,  (x,  y) 
in  image  spec*,  and  (  «,^)  in  the  exit  pupil.  The  cooruinates  in  image  space 
are  normalized  by  the  lateral  magnification  of  the  imaging  system.  This  is 
done  to  make  the  coordinates  of  a  given  object  point  equal  in  magnitude  to  those 
of  the  corresp  -noliig  image  point.  The  conventions  regarding  functional 
representation  are  : 

f(x)  a  f(x.y),  (5.1.1) 

^l'y2^  ’  (5.1.2) 

4Xj  ■  dXj  dyj  ,  (5.1.3) 

and 

X,  >  XjX.®  f  y^y^^  (5.1.4) 

Here  the  subscripts  1  and  2  denote  the  point,  or  P2,  whose  coordinates 
0  0 

are  used;  and  x  and  y  are  »init  vectors  in  the  direction  of  the  x  and  y  axis 
respectively.  The  same  conventions  of  course  apply  In  the  aoerture  and  object 
planes  ^8ee  Figure  3). 

We  shall  be  concerned  with  the  transmission  of  distributions  from  object 
apace  to  image  space  a!Kl  the  subscripts  0  and  1  will  denote  that  the  distribution 
is  an  object  01  image  re'ipectlvely,  i.e., 

1^(0  distribution  in  object  space  (5- 1-5) 

f.(x)  Corresponding  distribution  in 

image  space. 


FifSUas-  3 


Since  we  shall  require  the  Fourier  transform  of  space  functions,  we 
associate  with  each  space  coordinate  a  spatial  frequency  coordinate  using  the 
following  convention  :  is  the  spatial  frequency  associated  with  the 

Cartesian  C' ordinate  Xj  ,  and  pjy  is  associated  with  y^.  The  functional 
conventions  introduced  for  the  space  functions  will  of  course  also  be  used  for 
'he  spatial  frequency  functions;  i.e., 

t(tii)  ■  fhy)  ' 

0  0 

•til  "  ^  Miy  y  , 

and 

dHj  =  . 

tj 

y 


(5. I.e) 

(5.1.7) 

(5.1.8) 
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The  time  coordinate  is  denoted  by  t  and  the  time  delay  coordinate  by  t  . 

The  associated  temporal  frequency  is  denoted  by  v  . 

Contrary  to  the  convention  of  the  preceding  chapters,  the  mutual  coherence 
function  will  be  written  here  witli  its  full  argument  rather  than  with  the  sub¬ 
script  12. 

Associated  with  every  function  of  the  space  time  coordinates,  Ffepj^,  t  ) , 
will  be  three  other  functions;  namely,  its  "spatial"  Fourier  transform,  F(^,y2>T) ; 
its  "temporal"  Fourier  transform,  F(2Cj,X2,  i^)  ;  and  its  total  Fourier  transform 

o 

Fviijf  i.e.. 


F(Up  U2t'^)  = 


2^ri^■T 


(5.1. If/ 
(5.1.10) 


2fl(ypXj  +  P2-?2 

.  (5.1.11) 

Throughout  the  rest  of  the  chapter  all  integrals  will  be  written  with  a  single 
Integral  sign  without  limits.  The  order  of  the  integration  will  be  implied  by  the 
differentials. 

5.1.2  A  Generaliaed  Transfer  Function 

In  this  section  will  be  discussed  the  problem  of  determining  the  image  of 
an  bxtended  polychromatic  object  and  the  inverse  problem,  viz.,  that  of  determining 
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the  object  from  a  knowledge  of  the  image.  We  shall  show  that  the  solution  nf 
this  later  problem  is  fundamentally  impossible  with  systems  of  finite  aparture. 

The  object  is  considered  to  be  planar  and  is  specified  by  its  mutual  coherence 
function,  ro^^ri-2’^^  central  problem  is  to  formulate  the 

relation  between  Fq  ^  ^  in  terms  of  tlie 

aperture  illumination  function.  No  assumptions  concerning  the  spectral  width 
or  degree  of  coherence  will  be  made  in  this  section.  We  shall  show  that,  while 
the  evaluation  of  the  Integrals  might  prove  somewhat  formidable  in  certain 
practical  applications,  a  solution  in  closed  form  can  be  obtained  and  a  transfer 
function  defined  which  is  simply  related  to  the  pupil  function. 

From  the  linearity  of  Maxwell's  equation  i‘  follows  (see  Appendix  4)  that 
if  there  are  no  non-linear  devices  in  tlie  imaging  system  the  mutual  spectral 
density  will  be  propagated  through  the  system  in  accordance  with  two  linear 
differential  equations,  i.e., 

D,  [  1  ■  0  (..1,2)  ,  (5  1.12) 

where  D  i  a  linear  differential  operator  in  the  coordinates  of  .P  .  Solving 
“  *  A 

the  first  of  these  equations  we  obtain  the  mutual  spectral  density,  v)  , 

betw'een  the  oscillations  at  a  typical  object  point  and  those  at  a  typical 

image  point  .  Here  the  subscript  iO  denotes  that  the  function  depends 

on  a  point  in  the  object  space  and  a  point  in  the  image  space.  This  partial  eo- 

lut.on  may  be  obtained  by  using  only  the  linearity  of  (5.1.12)  as  follows  ;  Let 

the  contribution  to  the  "complex  disturbance"  at  g,  due  to  the  "disturbance" 

A 

from  an  element  d^,  of  the  source  around  be 

K(^j,Ij,  V)  d^  j  .  The  function  K(  v)  describes  the  optical  imaging 
system.  Then,  since  (5.1.12)  is  a  linear  differential  equation,  the  toUl  "disturbance" 
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at  is  given  by 


^q(5Sj»  4  2>  ^2’ ^^-1*  ip  (5.1.13) 

Next  v/e  repeat  the  argument  used  above  solving  this  time  the  second  of  the 
equations  (5.1.12).  For  the  sake  of  generality  we  assume  the  optical  sys¬ 
tem  to  be  characterized  by  a  second  function,  J(x2,  i2f^)  ■  Th®  relation 
between  J  and  K  will  be  determined  below.  Using  the  linearity  of  the 
remaining  equation  we  obtain  the  image  a,<? 


Substituting  from  (5.1.13)  into  (5.1.14)  we  obtain  finally 


(5.1.14) 


~  J(X2’ ^2’ ^^P^P  ^^^^1  *^^2 

. (5.1.15) 

Before  discussing  the  physical  significance  of  the  functions  J  and  K, 
we  shall  slvow  that  there  is  «  slnijj^e  relation  between  them.  To  this  end  we 
interchange  Uie  roles  of  and  Li  i.2  (5.1.15)  o'^tainlng 


(5.1.16) 


However,  by  Lemma  I  (Chapter  1) 
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and 


4  ^ 


Using  these  two  relations,  comparison  of  (5.1.16)  ana  (b  1.15)  shows  that 


. (5.1.17) 


From  (5.1.17)  it  follows  immediately  that 


(5.1.18) 


Using  this  result,  (5.1.18),  we  ma;/  rewrite  (5.1.15)  as 


A 

Kjii,^l,»^)K  (52'i.2’^^  ^ ^1  ^^2 

. (5.1.19) 


Equation  (5.^  19)  is  the  basic  relation  of  this  analysis.  It  expresses  the 
temporal  mutual  spectral  density  of  the  ima^e  in  terms  of  the  tempcral  mu- 
tr.al  spectral  density  of  the  object  and  a  function  K  which  characterizes  tne 
imaging  system. 

Before  continuing  the  development,  it  is  useful  to  consider  the  physical 
irti  rpretatlon  of  the  function  K(Xj.^,c).  The  demonstration  of  the  l.iterpre- 
tation  of  K(Xp  is  straightforward  but  requires  some  results  from  a  later 

section.  Therefore,  to  preserve  the  continuity  of  the  present  discussion,  this  demon¬ 
stration  is  given  In  Appendix  3  and  we  simply  state  here  that  the  function  K 
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represents  the  complex  amplitude  at  Xp  due  to  a  monochromatic  point  source 
at  of  frequency,  v. 

Until  new  we  have  not  specified  the  position  of  the  image  plane.  Since  there 
are  by  -ieflnidon  no  imaging  elements  between  the  exit  pupil  and  the  image  plane 
it  is  clear  that  equation  {5.l.x9)  is  valid  throughout  this  entire  region.  However, 
since  the  form  of  K  in  the  exit  pupil  varies  significantly  from  its  form  in 
Gaussian  image  plane,  and  since  its  behaviour  in  each  of  these  planes  is  of 
particular  physical  importance,  it  will  prove  convenient  to  designate  by  two 
different  symbols  Uie  form  of  K(Xp  on  these  two  surfaces.  Accordingly, 

we  shall  denote  by  A(  i u)  the  form  of  K  in  the  plana  of  the  exit  pupil.  Thus 
A(i,  ii,  I')  is  complex  disturbance  ai  a  point  in  exit  pupil  due  to  the  mono¬ 
chromatic  point  source  of  frequency  at  a  point  in  the  object  plane.  We  reia. 
the  symbol  Ki^^,  v)  to  denote  the  function  in  the  image  plane.  Using  this  c  on¬ 
vention,  K{x,  4,  v)  may  be  thought  of  as  tl  distribution  in  the  image  plane  due 
to  a  monochromatic  distribution  of  complex  amplitude  in  the  plane  of  the  exU 
pupil  rememberl!^  of  course  that  the  distribution  in  the  exit  pupil  i.s  determined 
by  the  object.  The  relation  between  those  two  functions  may  then  be  expreseedas 

K(x,i,  -  J Mi.  ^  (5.1.20) 

Here  a  is  a  point  in  the  aperture  of  the  system;  G  i-  a  Green's  function  satis¬ 
fying  the  Helmholtz  equation  and  vanishing  over  the  plane  of  the  exit  pupil. 

Under  the  conditions  characterizing  most  imaging  systems  (5.1.20)  takes 
a  particularly  simple  :..m;  but  before  discussing  this  point,  we  shall  continue 
the  development  up  to  the  Introduction  of  the  general  transfer  function.  This  is 
done  to  avoid  the  erroneous  impression  that  the  transfer  function  analysis  ^ 

y 

I 

I 
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involves  the  Fraunhofer  (or  far  field)  approximations,  as  is  sometimes  be¬ 
lieved  to  oe  tile  cace. 

Equation  (5.1.19)  assumes  a  convenient  and  useful  form  if  the  system 
under  consideration  is  "spatially  statio.iary",  i.e.,  if  the  function  is 

a  lui  ",tlon  of  the  difference  of  the  spatial  coordinate j, 

K(x,i,  v)  a  K(i-5,  ^)  ,  (5.1.21) 

This  condition  is  satisfied  by  scanning  systems,*  which  includes  most  an¬ 
tenna  systems  and  also  by  many  Important  (vislolc)  optical  systems.  For 
systems  which  do  not  scan  but  form  the  ontirt  spatial  image  slmulianeously,  K 
will  not  in  general  be  a  function  of  the  difference  (i  -  i)  only.  Howe  a r,  for 
most  optical  systems  tlie  form  of  the  diffraction  pattern  varies  slowly  across  the 
image  plane.  Hence,  the  image  space  may  be  divided  into  "isoplanatic"  areas 
over  which  K  may  be  aasumed  to  be  a  function  of  (^  -  5)  to  any  desired 
accuracy.  This  posslb-Uity  Is  discussed  at  length  by  Fellgett  and  Linfoot  (1955) 
and  by  Duraontet  (1955). 

Tliroughout  the  rest  of  thio  discussion  wf  shall  0  ly  be  concerned  wi?h  systems 
foi  which  the  condition  (6.i.2l)  is  satisfied. 


By  this  we  inean  Sj stems  in  which  the  image  forming  device  (e.g.  antenna) 
scans  while  the  position  of  the  detector  (e.g.  feed)  relative  to  the  aperture 
remains  fixed,  I.e.,  the  antenna  and  feen  move  together.  Tlie  above  con¬ 
sideration  Is  not  valid  for  a  system  in  whlc.h  the  Image  forming  device  re¬ 
mains  fixed  while  the  detector  scans  the  aerial  Image.  A  camera  with  a  focal 
plane  shutter  is  a  simple  example  of  such  a  system;  the  lens  creates  the 
entire  im.ige  simultaneously  and  the  shutter  slit  then  scans  the  image. 


Using  the  stationarity  corsdition,  (5.1.19)  may  be  rewritten  as 


K(ii  -  Sp  t-)  K*  (ij  -  1>  ^2’  *^il  ^^^2 

...  .  (5.1.22) 

Taking  rr/«  "snace-type"  Fourier  transform  of  both  sides  of  (5.1.22)  and 
using  the  convoluiior;  theorem  yields 

0  0 

=  K(iipi^)  X  ■  (5.1.23) 

At  'hia  pc^ni  w  ,  introduce  the  transfer  function,  L  1^2' defined  by 

X  K(iip^)  ("Jig,  '  p)  -  (5.1.24) 

L  (£.1.2)  becomes  simply 
0  0 

(5.1.25; 

Equation  (5.1.2)  may  be  resardeu  as  the  basic  equation  in  tJ»e  frequency 
cfcujain  analysis  of  lms;£ing  systems,  it  is  clear  from  the  foregoing  dlacusslon 
that  the  transfer  fsinrtion  analysis  is  rigorously  applicable  to  any  "spatially 
stationary"  system.  That  is  no  approximation  need  be  made  ccncerni.ng  the  re¬ 
lation  between  the  aperture  illumination  function,  A  ,  and  the  diffraction  pattern, 
K-  Howover,  in  tnost  optical  applications  and  in  fact  in  njost  antenna  applications, 
the  diffraction  patieni  is  characterized  by  ^he  Fraunhofer  approximations  (or 
thtr  foimaliy  equivalent  far  field  approximations);  and  since  under  these  condi¬ 
tions  (5,1.20)  assumes  a  particularly  simple  form,  we  shall  introduce  these 
approximations  at  this  point  and  retain  them  throughout  the  subsequent  sections. 
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Under  tiie  usual  approximations  whitj'i  characterize  i  .aunhofer  diffraction, 
(5.1.20)  reduces  to 

f  -2jr  ivvcT  •  i  ) 

(5.1.26) 

J 

where  X  is  the  wave  length  of  tiie  spectral  component  belonging  ro  frequency 
i>  and  R  is  the  radius  of  the  Gaussian  reference  sphere.  However,  K(x,^,i^) 
may  be  expressed  as 


«  Jfe,  a,  Sii.. 


K(x,i.^)  -  I  K(x  (5.1.27) 

By  comparing  (5.i.26)and  (5.1,27)  we  Identify  the  spatial  frequency,  ji,  as  the 
reduced  aperture  coordinate,  i.e.. 


Further  we  note  the  important  relation 


K(x,iii')  =  K(j,  ^,c)  =  A(x,ii^)  . 

Our  aubsfl'.,  lenl  analysis  will  deal  soiely  with  K  and  Its  transform  I<  From 
(5.1.24),  (5.1.26)  and  (5.1.27)  It  is  clear  that  under  the  conditions  stated  the 

transfer  function,  L(i)pii2«*'^'  *  systen.  utilizing  partially  coherent  poly¬ 

chromatic  illumination  is  the  product  of  the  frequency  dependent  aperture  illumina¬ 
tion  function  considered  as  a  fvuictlon  of  spatial  frequency  and  evaluated  at 
t^l  =  flIj/'XH  with  Its  complex  conjigatc  evaluated  at  ^2  =£f2A^- 

It  is  clear  from  the  above  considerations  that  the  inverse  pr-oblem,  namely 
that  of  determining  the  object  from  a  knowledge  of  the  Image,  canr'ot  be  solved 
if  the  Imaging  system  has  a  finite  aperture.  This  roncUisto.n  can  t>e  under^too-u 
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V’V  formally  inve”tmg  (5.1.25).  We  then  obtain 

ro^yp  1!:2>  X  0  .  (-01.23) 

X  (l^pM2.'^) 


However,  from  (5.1.27)  and  (5.1.24)  it  followc  that  if  the  aperture  iS  finite 
^  is  identically  zero  beyond  some  maximum  frequency  j  u  j  max,  and  hence, 
(5.1.28)  is  indeterminate.  Thus  the  inverse  problem  is  soluOle  only  up  to  an 
arbitrary  function,  f  , 


0 

rj(l^pii2> 


(5.1  29) 


where  f  is  any  function  of  frt  juencies  greater  than  j  n  j  . 

Equation  (5.1.25)  gives  the  relation  between  the  total  spec ‘ral  densities 
of  tile  obj^'ct  and  image  and  is  thus  the  solution  sought  in  this  section,  ^lowever, 
in  many  applications  one  is  interested  in  a  much  less  general  solution,  namely, 
the  intensity  distribution  in  the  image.  Obtaining  the  intensity  distribution  in 
gen  ral,  from  (5.1.22),  is  somewhat  involved  and  not  very  helpful.  However, 
in  the  limiting  cases  of  coherently  and  incoherently  illuminated  objects  the 
problem  is  tractable,  and  in  the  subsequent  sections  we  shall  discuss  these 
limits  in  detail  for  botli  polychromatic  and  quasi-monochromatic  illumination. 

5.2  The  Limiting  Forms  of  the  Transfer  Function 

Following  the  pattern  established  in  die  earlier  chapters  we  first  examine 
the  limiting  forms  for  polychromatic  light  and  in  a  later  section  examine  the  ex¬ 
tremes  under  the  quasi-monochromatic  approximation. 

We  shall  show  in  this  section  tliat  the  ♦ransfer  function  for  coherent  objects 
is  the  frequency  dependent  aperture  Illumination  function  K  -  A  (x,  i^) 


* 

e  . 


.4  ^ 


evaluated  at  i^q  and  considered  as  a  function  of  spatial  frequency.  (Here  Vq 
is  the  frequency  of  the  illun'ination.  See  Chapter  3^ 

TTie  analysis  of  the  systems  imagii  g  incoherent  objects  is  complicated 

by  the  fact  that  as  explained  in  Chapter  4  the  image  is  partially  coherent. 

Thus,  if  one  seeks  a  complete  solution  (i.e.  the  mutual  coherence  of  the  image), 

the  ti*an8fer  function  must  operate  on  the  spectral  density  of  the  object,  l{^v) 

(a  function  of  one  point  only)  to  produce  the  mutual  spectral  d  nsity  of  the  image, 
0 

function  of  two  points).  This  consideration,  overlooked  or 
omitted  in  the  literatvme  on  the  imaging  problem,  is  important  in  the  treatment 
of  cascaded  systems.  Ihe  required  transfer  function  will  be  shewn  to  be  the 
function  L  introduced  in  the  previous  section. 

If,  on  the  other  hand,  one  requires  only  the  Intensity  distribution  in  the 
image,  l(j|,0),  the  entire  analyola  may  be  performed  wiLh  functions  of  one 
point  only,  the  i^ctral  donoltlea  of  the  object  and  tmage.  For  this  problem, 
the  ‘ransfer  funcU^^n,  M  ,  will  be  ehown  to  be  the  convolution  of  the  frequency 
depeiKient  aperture  illumination  with  its  complex  conjugate. 

5.2,1  The  oherent  Limit 

It  was  shown  In  Chapter  ^  in  a  coherent  field  the  mutual  coherence 
function  is  of  tht;:  form 


.  -2»iV«T 

Todi-la-')  ' 


(5.2.1) 


With  a  mutual  spectral  density 


i5(5i.S2’‘^  ■  W'o  (52>  »  f-  -  *'o» 


(6,2.2) 


Sutotituting  from  (5.2.2)  into  the  general  solution,  (5.1.19),  and  tpking  the 
temporal  Fourier  transform  on  both  sides  we  obtain 


where 


*  -21ri^»»T 

UjCSi)  ^i 


(|)  K  (  I  -  x)d{ 


(5.2.3) 


(5.2.4) 


From  (5.2.3)  and  the  theorems  of  Chapter  3  it  follows  immediately  that; 
Theorem  XIV  :  The  image  of  a  coherent  object  is  coherent. 

Taking  the  spatial  Fourier  transform  of  both  sides  of  (5.2.4)  and  using 
the  convolution  theorem  we  obtain 


K  (iij,V(j)  .  (5.i.5) 

The  appropriate  transfer  function  is  the  frequency  dependent  aperture  Illumina¬ 
tion  function  evaluated  at  Pq  and  considered  as  a  function  of  spatial  frequency, 

“  A(^jj,Vq).  The  coherent  image  is  completely  determined  by 
(5.2.5),  and  ihe  intensity  distribution  Is  obtained  as  a  special  case  of  (5.2.3)  by 
setting  5^  =  X2>  t  =  0. 

5.2.2  The  Incoherent  Limit 

An  incoherent  object  is  described  by  a  mutual  coherence  function  of  tJje 


form  (see  Chaptes  3) 


A 


n 

%liere  IqCI^,  i)  is  the  bsH  coherence  function  sit  defined  by  the  relation 

I(i,T)  =<V(|,t  +  T)  V*(i,t)>  .  (5.2,7) 

The  mutual  i^ectral  density  is,  therefore, 

=  ’(ia-*')  >  <ia  -  j'')  ■ 


Substituting  from  (5.2.8)  into  the  general  aolut^on,  (5.1.19),  w«  obtain  the 


1f 


K(|i  K  (|2  -*2'*'^  ^^2“  **^1  ^2 

. (5.2.9) 


We  may  now  integrate  over  d^^  obtain 


n(Sl.*«.»')  =  K*(^j-g2,»')  l(ii,>')dii  ;  (6.2.10) 

and  taking  the  spatial  Fourier  transform  of  both  side:  and  using  again  the  convc' 
lutlon  theorem  yields 


0 

“  ‘^1^1  M2»  ^2'*'^*  (5.2.11) 


where  Is  the  gene raibed  transfer  fuse tlon  defined  in  Section  5. 1 

^  (iLvUtt  *')  ■  K(+iij,  +1')  .  .  (5.2.12) 


f 


It  is  convenient  to  rewrite  (5.2. 11)  in  the  form 


U  A 

•  (5.2.13) 

Equation  (5.2.  i3)  eitpresses  the  fact  that  the  energy  contained  in  the  object  »' 
at  fi  is  distributed  among  aU  pairs  of  frequency,  and  j£  -  ki  the  image. 

While  equation  (5.2.13)  is  the  complete  solution  (for  an  incoherent  object) 
for  the  total  mutual  spectral  density  <n  the  image  and  hence  for  the  mutual 
coherence  function,  one  is  often  interested  in  the  more  restrictive  solution,  the 
intensity  in  the  image.  This  is  obtained  at  once  by  settii^  Ij  =  JSj  in  (6,2.10), 
which  gives 

-  J |K(^j  .ij,i^)|  ^I(^j,v)d^j  . 

.  (6.2.14) 

Taking  the  spetial  Fourier  tranalurm  of  both  sides  of  (S.2.14)  we  obtain 

->l(ii,v)  ,  (6.2.16) 

where 

'>n  (if,!')  =J\k($,v)  l^e  .  (5.2.18) 

From  thT  interpretation  of  K  given  in  Appendix  4  it  is  clear  that  |  K  |  *  is 
the  frequency  impendent  intensity  diffraction  pattern  oi  the  lmagii«  syatem.  The 
transfer  function  for  determining  the  total  spectral  density  of  the  image  of  an  in¬ 
coherent  object  is  thus  the  transform  of  the  "intensity  diflractkin  pattern"  of  Ufe 
system.  The  formula  (5.2.16)  is  simplified  further  by  again  using  the  convolidlon 


ti'.eorem  find  (5.1.2)  which  gives 


J 


'hL(^v)  »  /  K(g-|i,v)  K*(g,i/)dc; 


(5.2.17) 


Hie  tpfttial  epectr&l  density  may  nov/  be  obtained  by  taking  the  Fourier  trans¬ 
form  of  both  sides  of  (5.1.15)  and  evaluating  at  r  =  0  ;  thus 


(5.2.18) 


The  formula  (2.2. 18)  eiqiresses  the  tect  that  each  temporal  spectral  component 
contributes  separately  and  independently  to  the  energy  in  the  spntial  frequency 
component  |i .  The  intensity  distribution  in  the  image  is  tlien  gb'en  by  the  spatial 
transform  of  (5.2.18). 

5.3  Imaging  with  Quasi- Itonochroma tic  Light 

While  many  Imaging  systems  of  practical  and  theoretical  Interest  deal 
with  polychromatic  light,  only  tiae  problem  of  imaging  with  quasi-monochromatlc 
or  monodiromatic  illumination  has  been  extensively  discusjied  in  the  literature. 
This  omission  of  the  more  general  problem  is  sssily  understood  since  without 
a  rigorous  and  general  formulation  of  coherence  theory  the  noathematlcal  analyaig 
Is  prohibitive;  and  the  earlier  formulations  of  this  theory  diacusiied  In  Chiqtter  1, 
are  not  well  suited  for  extension  to  this  general  caee.  Hiei-efore,  the  transfer 
function  analysis  as  found  in  the  svsilslde  literature  is  applicable  only  to  quar i- 
monochromatic  light,  and  in  order  to  compare  the  results  given  hero  with  those 
of  earlier  writers,  ws  examine  in  this  section  the  form  ol  tue  transler  functions 
under  the  quasi- monochromstic  aiguoximstion. 

A  psrtlalty  coherent  qussl-monochromatic  object  will  be  described  by  a 
mutual  ooherenee  fuMtion  of  die  form  (see  Chapter  3) 
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-2t  iVT 


ro(i,,i2,r)  ^  r5(i,,i2.»)®  (lTi<<^):,  (5.3.1) 


its  temporal  transform  is 


Substituting  from  (5.3.2)  into  the  goneral  solution  (5.1.19)  and  taking  the 
inverse  Fourier  transform,  we  obtain 

A  ••27riPT 

. (5.3.3) 

Equation  (5.3.3)  provides  the  starting  point  for  the  analysis  of  partially  co¬ 
herent  images.  Beginning  from  (5.3  3)  the  entire  atialysis  of  the  two  pre¬ 
ceding  sections  may  be  taken  over  mutatis  mutandis  for  the  quasi-mnnochromatic 
imag'ng  problems  considered  here.  DenQi...ig  by  the  suffix  q  that  the  functions 
are  applicable  to  quasl-monochromatic  light,  the  various  transfer  functions 
and  frequency  domain  image  equations  are  : 

Coherent  Object 

transfer  function  'K  a  ~ 
imaging  equation  Uj(ii,  P)  -^q(ji,P)  Ug(ii,P) 

Panlally  Coherent  Object 

transfer  function  ^  q  =  K* 

imaging  equation  H (ifj, P2,»>)  =  ^q(ypy2.P) 

I 


/ 


K  d|,d£2 


4^ 

Incoherent  Object 

transfer  function  ajq  '  K( 

imaging  equation 

transfer  faction  =  /  KKa~y,v)  }i*(q, 

imaging  equation  5)  v) .  IqU  P) 

The  transfer  functions  for  quasi- mo t^ochromatic  light  are  thus  seen  to 
be  simply  the  generalized  transfer  functions  evaluated  at  the  mean  frequency. 
This  result  is  to  hav^.  been  e^qpected  since  mathematically  the  quasi- monochro¬ 
matic  approximation  ie  characterized  by  an  approximately  monochromatic  mutual 
coherence  function.  However,  it  should  be  emphasized  tliat  while  the  ti-ansfer 
functk  ns  obtained  in  this  section  may  be  formally  obtained  by  taking  a  single 
spectral  component  of  the  general  solutions  the  Inverse  procedure  (integrating 
the  quasi-monochromatic  solution  over  frequency  to  obtain  the  general  solutions) 
is  not  justifiable.  This  conclusion  is  evident  from  the  fact  that  the  solutions 
obtained  in  this  section  are  only  approximate  {  |  rj  <  <  and  accordingly 
the  transfer  functions  depend  on  the  mean  frequency,  v  ,  not  on  an  isolated 
frequency,  .  It  Is  thtis  logically  impossible  to  obtain  the  general  so¬ 
lutions  from  the  quasi-monochromatlc  solutions. 

Hie  transfer  functions  defined  in  this  chapter  are  summarized  in  Table  2. 

N.B.  Two  imaging  equations  are  required  for  systems  Involving  in¬ 
coherent  objects.  is  u>ed  to  determine  spatial  spectral  density 

in  the  image,  while  th#  more  general  problem  of  determining  the  mutual 
spectral  density  involves  the  transfer  function  (jaj,  2  »^)  •  However, 

both  of  these  functions  operate  on  the  spatial  spectral  density  of  the  object. 
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APPENDIX  1 


iiSLATION  BETWEEN  VON  LAUE'S  MEASURE  OF 
COHERENCE  AND  THE  COMPLEX  DEGREE  OF  COHERENCE 


TTie  measure  of  coherence,  ,  introduced  by  von  Laue  (1907)  is 
defined  by 

y  --  - — M — — — n - .  (A.  l.i) 


Here  fj(t)  and  £2(1)  are  real  funciknn  of  time  which  describe  the  (scalar) 
optical  distKvbonce  at  the  two  pjinta  Pj  and  Pj » respectively.  The  function ' 
^j(t)  (i  =  1,2)  ha/e,  therefore,  the  Fourier  representation 


K 

tl(t)  =  J  Fj  (V)  coe[^j(r)  -  iwvt  ]d»/  , 


(A.  1.2) 


and  tue  functions  gj(t)  are  defined  as 


!,«)  =  /  F, 


(;^)8ln  {  ^(t')  -  2t>'t]dt'  . 


(A.  1.3) 


This  formulation  strongly  suggests  the  introduction  of  the  analytic  signal; 
and  accoruiogly,  we  may  rewrite  (A.I..)  as 


<Vi''  (tiXv/  (l)> 

where  the  functions  V-*’  and  are  Identlllwl  as  t  and  g,  re^ctiveiy. 


(A.  1.4) 


i 
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In  order  tc  obtain  the  relation  between  and  we  have  but 

to  recall  the  definition  of  y22(‘^)  •  Equation  (1.2.10)  of  Cluipter  1  may 
be  formally  e^qumded  to  give 


y,2(T)  -- 


(A  +  iB) 


j(t)Vj*(t+T)>  <V2(t)V2"{t+T)S 


where 


A  -  ^  T)y  +  <0'At)V2^(t  >  T)} 


B 


[<Vj‘(l)Vj’'(t+lj>  - 


(A.1.S) 


but  by  the  theorema  eatabliiihed  in  Chapter  2 


and 


<V2*‘(t)  =  <v/(t)  V2‘(t4T)^  , 

<^v/(t)  =  -<(Vj*’(t)Vi‘(t+T^, 


whence  Equatton  (A.  1.5)  becomes 

a  [<v/(t)  V3'’(t+T)>4  icVj^t)  v,*‘{t^Tj>] 

fj  j(  >1  ”  ? 


/  <V,(t)  v,*(i.  j>  <(v,(t)  v,Vti> 

Multiplying  both  BkiSS  01  ^  and  aoUj^f  from  (1.>.14)  that 


<^Vj(g  Vj*a+T))>  -  2  <7/(1)  Vj^(t+T)^ 


I 


36 


we  obtain 


IVliCr) 


<(v^^  (t)  V2^  (t + 7i>  2  +  <y/  it)  V2"  (t + Tp 
<(vi'’  (t)  {t+  T)^  <;^  Vg*’  (t)  Vg**  (t  ^ 


(A. 1.7) 


or  at  T  =  0 


(t)  v/  t  ^Vj'  (1)  v/  (tj>  ^ 
<y/^  (i)><v/‘  (tj> 


(A.  1.8) 


Compariaonol  (A- 1.8)  and  (A.  1.4),  yields  the  relation  between  von  Laue' a 
measure  of  coherence  and  Wolfs  complex  degree  of  coherence  function  as 


(A  1.9) 


APPEJ®IX  2 

THE  QUASI- MONOCHROMATIC  APPROXIM.\TION 
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We  may  define  quasi-monochromatic  illumination  by  the  following 

A 

property  :  the  mutual  spectral  density,  appreciably  different 

from  zero  only  for  those  sp^ct-^al  conjponents,  which  satisfy  the  inequality 

[f  -  P  I  «  Ai^  , 

where  F  is  the  mean  frequency  and  A  is  the  spectral  width  of  the  lignt. 
Physically  this  cor>dition  implies  that  most  of  the  energy  in  the  field  is  con¬ 
tained  in  the  spectral  reg  on  r'-Ai'<i'vF+Ai'' 

Tc  obtain  the  form  of  the  mutual  coherence  function  for  quasi-mono- 
chromatic  fields  we  first  recall  that  t>o  expressed  in 

the  form 


A  -2Trli^T 

n„(i^)e  d^-  . 


(A  2.1) 


We  may  now  factor  out  the  mean  freqiency  teim  of  the  integrand  in  (A.2.1) 


and  obtain 


-2  iriPr  r  -2  ffi(p-  P)  r 

®  /  ris(‘'')e  dr  .  (A2.2) 


If  we  limit  our  attention  to  sufficiently  small 


more  precisely  if 


■j- 1  «  -i-  the  frequency  dependent  factor  of  the  exponent  of  tl'e  integrarvl 

!  A  L/  ^ 


satisfies  the  inequality 


(r-p)  T  « 


for  all  values  of  v  for  whlr:h 


,(v)  I  is  significant.  Thus,  the 


variations  in  e  ^  may  to  3  good  approximation  be  ignored,  and 

(A  ''.2)  may  be  rewritten  as 


n  ,'.A  - 

i  12* 


-2  nivT 


r,2(v)d^  (!t  «  i)  .  (A.2.3) 


The  integral  in  (A.2.3)  may  be  formally  evaluated  to  give 


rn,(.o 


r,9.(o)  . 


and  we  have  finahly 


^2^^  -ri2(0)e 


(A.2.4) 
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APPENDIX  3 

THE  GENERAL  UNIMODULAR  ANALYTIC  SIGNAL 
AS  AN  AUTOCORRELATION  FUNCTION 


We  will  show  in  this  appendix  that  the  general  unimodular  analytic  signal 
can  be  interpreted  as  an  autocorrelation  function  only  in  the  degenerate  case 
that  it  can  be  written  as 

e  =  e 

The  constants  a^^*  determine  the  position  of  the  poles  of  the  meromorphic 


function 


Aj  e 


2^ 


'  00 

- 

TT  !n 

f  O 

n=l  a^ 

a  *-  z 

In  1 

i03  -  2irv  qZ) 


(A.3.1) 


In  (A.3.1)  the  imaginary  part  of  has  the  same  sign  as  i.e., 


l{aj^|<  0,  (A.3.2) 

and  ^  and  i/q  are  real  constants. 

The  most  general  unimodular  analytic  signal  la  obtained  by  setting 
z  =  T  in  (A.3.1)  (see  Chapter  2).  The  requirement  that  (A.3. 1)  repre- 
pei't  an  autocorrelation  function  Is  equivalent  (see  Chapter  3)  to 


il 


90 


00 

f  fiiW 


sini/ri/TdT 


3  0  , 


(A.  3. 3) 


where 


2 

fll(T)  -  e 


(A.  3.4) 


We  shall  show  here  that  (A. 3. 3)  can  be  satisfied  if  and  only  if  the 
meromorphic  function,  (A.  3.1),  has  ro  poles. 

The  integral  in  (A,  3. 3)  is  conveniently  evaluated  in  the  complex  plane. 
Equation  (A.  3.1)  can  oe  v/ritten  as 


where 


R  +  F_]  =  0  , 

OO 

,  f  +2jr  ii^^z 


(A.  3. 5) 


(A.3.G) 


In  (A.  3,5)  and  (A.  3. 6)  the  identity, 


sin  2-n 


?.n  iVpZ  -2ff 

e  -  e 


(A.  3. 7) 


was  used.  The  function  in  (A. 3. 6)  Is  given  by  an  integral  aloi^  the 

.’Gal  axis,  i.e., 


T? 

-  Tr 


n:--l 


♦  / 


z\]  z] 


(A  3.8) 


i 


Hie  integral  in  (A.  3.8)  can  be  evaluated  by  contour  integralion  closing 
the  contour  at  infinity  above  or  below  the  real  axis  accordingly  as  Vq  +  v 
is  less  than  or  greater  tlian  zero.  Thus  by  Cauchy's  residue  tlieorem 
0  u<- 


E  I  i 

Tiff  m=l  n=l  a 
n^m 


I  a  -  a  t 

n  n  m 

I  \  a  -  a  *  I 
n  \  n  ml 


*1  ip 
*16  6 


'  “  >’>-  ■'oj. 

. (A.  3.9) 


since  by  (A.3.2)  the  poles  are  all  7n  tiie  lower  half  plane  Similarly, 


0 

oo  oo 


<  Vi 


.  }a"  E  TT  \  hn  .2ffi(vp~v)a^  \ 

-\Thf  n.=l  n=l  ■£“  r  ^  r  ^ 

n?-'m  1  n  m  / 


(A.  3. 10) 


Tne  product 


I  T  a  /  -  z  \ 

1  !  _n  _n  _ 

n-i  -  z 

n  \  n 


(A.  3. 11) 


convei-ges  for  aii  z  (c.t.  Chapter  2).  li  tl.erefore  converges  for  z  - 
and  we  may  write 


'a  - 

a 

n 

in 

a  ♦  - 

a  * 

n 

ni 

=  ^  IC^ 

m  xn 


(A.  3. 12) 
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Using  (A.  3. 12)  we  may  write 


F+  + 


^  (B^  ^‘^m^  t  ^  i  sin  2n(i^  -1-1/  Q)bjJ  e 


27r(v+i/Q)C 


m 


m=l 


+  2  (B„  + 


m=l 


m  m' 


cos  27r(i'Q-i/)b  -h  i  sin  27r  (i/q  +p)  b, 


1  27r{i/Q-J/)C 


m 


m 


. (A.  3. 12) 

In  (A.3.13)  and  are  the  r-^al  and  imaginary  parts  respectively  of  a^; 
and  *B  greater  than  zero  (see  Chapter  2)  The  physical  requirement, 

(A. 3.1),  Is  that  +  F_  vanish  identically  tor  all  v.  Equation  (A.3.13) 
can  be  written  in  a  more  convenient  form  as 


)  cos[<^j^+2]r(;/Q-i/)bJ«  0  i/>i/o 


where 


R 


m 


®in^  ^  ^  1b~  ^ 

m 


(A.  3. 14) 


(A.  3. 15) 


That  this  equation  cannot  be  satisfied  for  R^^  0  is  clear  from  the  asymptotic 

behaviour  of  the  sums.  The  second  sum  vanishes  for  large  v  while  the  first 
sum  diverges  term  by  term  as  e*^  .  Thus  we  must  have 


■  0 

m 


(A.  3. 16) 


However,  since  In  the  solution  of  the  integral  equation  (3.7)  poles  at  the 
origin  were  excluded  and  the  product  was  taken  only  over  the  poles  occurring 
for  finite  z,  the  cannot  vanish  if  there  are  poles  or  zeros  in  the 
meromorphic  function.  Consequently  there  can  be  no  poles  for  finite  z  if 
(A.  3.1)  is  satisfied. 


Thus 
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EXISTENCE  OF  A  LINEAR  DIFFERENTIAL  EQUATION 
FOR  THE  MUTUAL  SPECTRAL  DENSITY 

From  Maxwell's  equations  it  follows  that  each  spectral  component  of 
the  optical  disturbance,  (t)  ,  satisfies  a  linear  differential  equation,  in 
particular  the  Helmholtz  equation  in  each  medium  comprising  the  imaging 
system.  Hence  we  may  write 


D  [  ^  (y)  ]  -  D[llm  V^(T,i^)]  -  0  ,  (A.4.1) 

T-^ 

where  V  (T,  i')  is  the  spectrum  of  the  truncated  function  Introduced  in 
Chapter  2  and  D  is  a  liuear  differential  operator.  Using  {A.4. 1)  and  the  con¬ 
volution  theorem  of  Fourier  analysis  we  may  write 

«0 

11m  D  [  (T,  V)  ]  =  lim  r D  [  V  (u)  ]  dv  =  0, 

T-«  T^  (  2P(^^-T) 


or  simply 

D[V‘'(T,»^)  ]  -  0.  (A.4.2) 

Thus  each  spectral  component  of  the  truncated  function  satisfies  the  saiue 
differential  equation.  That  the  function  likewise  introduced  in 

Chapter  2,  also  satisfies  the  same  equation  can  be  seen  as  follows.  By 
definition 

«o  « 

J  (*■) '  x*'  ■  T  J 

-to  -to  -to 


2  vlvt 


V  tT,v)e 


dwlt' 


1 


f 


i 


A 


> 

'I'' 
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Operating  on  both  sides  with  D  and  using  (A.4.2)  we  find 

D[tV^(v)  ]  =  0.  (A.4.3) 

From  the  linearity  of  D  and  equations  (A.4.2)  and  (A.4.3)  it  follows  that 

D[  V(T|i;)  ]  =  0,  (A.4.4) 

where  V(T|  v)  -  (T,  i*)  +  i  is  the  spectrum  of  the  analytic 

signal  associated  with  the  truncated  function  V^(T,y).  Equation  (A.4.4) 
expresses  a  particularly  convenient  property  of  the  analytic  signals,  vis., 
that  an  analytic  s^nal  sati.sfies  the  same  linear  differential  equation  as  the_ 
real  function  with  which  it  is  associated. 


We  recall  from  Chapter  2  that 


Vj(Tli/)V2  (T 


2T 


and  operating  on  both  sides  with  D^,  (s  =  1,2),  where  the  subscripts  denote 
that  the  operation  Is  in  the  coordinatea  of  P  we  obtain  finally 

D 


A 

f  rj2  -  0  (  8  =  1,2) 


(A.4.5) 


Thus  the  mutual  spectra’  density  will  be  propagated  through  the  system  by  the 
same  equations  which  govern  the  optical  disturbance  Itself. 
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APPENDIX  5 

INTERPRETATION  OF  THE  RESPONSE  FUNCTION  K(s,  i,v) 

It  was  shown  in  Chapter  4  that  a  coherent  (scalar)  field  may  be  com¬ 
pletely  described  by  a  single  conaplex  wave  function  U.  In  Chapter  5  it  Is 
shown  that  this  complex  function  is  of  central  importance  in  treating  the  image 
of  coherent  objects.  To  determine  the  physical  significance  of  we 

consider  a  coherent  point  object  as  the  input  to  the  imaging  system.  For  this 
object  the  function  U  is  of  the  form 


(A.5.1) 

Substituting  from  (A.  5.1) 

into  (5.2.4)  and  integrating  over  d 

we  obtain 

K(^q-^,»')  6(^'-^'o)  , 

(A.  5. 2) 

or 

(A.  5. 5; 

'Hie  physical  Interpretation  of  K(  is  clear  from  (A  f  .S);  that  is 

w)  is  the  complex  amplitude  in  a  diffraction  pattern  due  to  a  point  source 
of  frequency  v  located  at 
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